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Conversion  Factors 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as 
follows: 


Multiply 

By 

To  Obtain 

acres 

4,046.873 

square  meters 

cubic  feet 

0.02831685 

cubic  meters 

cubic  inches 

0.00001638706 

cubic  meters 

degrees  (angle) 

0.01745329 

radians 

degrees  Fahrenheit 

(5/9)  x  (°F  -  32) 

degrees  Celsius 

degrees  Fahrenheit 

(5/9)  x(°F- 32) +  273.15. 

kelvins 

feet 

0.3048 

meters 

gallons  (U.S.  liquid) 

0.003785412 

cubic  meters 

horsepower  (550  ft-lb  force  per  second) 

745.6999 

watts 

inches 

0.0254 

meters 

kips  per  square  foot 

47.88026 

kilopascals 

kips  per  square  inch 

6.894757 

megapascals 

miles  (U.S.  statute) 

1.609347 

kilometers 

pounds  (force) 

4.448222 

newtons 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass) 

0.4535924 

kilograms 

square  feet 

0.09290304 

square  meters 

square  miles 

2,589,998 

square  meters 

tons  (force) 

8,896.443 

newtons 

tons  (2,000  pounds,  mass) 

907.1847 

kilograms 

yards 

0.9144 

meters 

Systeme  International  d'Unites  (“International  System  of  Measurement”),  commonly  known  as  the  “metric  system. 
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1  Introduction 

Background 

In  the  Department  of  Defense,  encroachment  refers  to  any  activity  that  limits 
the  ability  to  continue  training  and  testing  at  military  installations.  Many  en¬ 
croachment  challenges  are  related  to  the  development  of  residences  near  mili¬ 
tary  installations — a  process  called  urbanization.  Urbanization  patterns  develop 
in  direct  response  to  past  regional  investments  in  infrastructure  (e.g.,  roads, 
sewer,  electricity,  and  gas)  and  regional  policies  (e.g.,  zoning  and  tax  incentives). 
The  free-market  system  then  allows  people  to  compete  for  land  and  to  develop  it 
within  the  opportunities  afforded  by  the  local  economy  guided  by  the  regional 
investments  and  policies.  Therefore: 

•  The  current  opportunity  to  train  and  test  is  the  result  of  past  regional  plan¬ 
ning. 

•  Future  opportunities  to  train  and  test  will  be  the  result  of  current  and  future 
investments  and  policies. 

It  is  important  that  proposed  investments  and  policies  be  evaluated  with  respect 
to  the  impact  on  future  urban  patterns  and  those  patterns  be  evaluated  with  re¬ 
spect  to  changes  in  opportunities  to  conduct  military  training  and  testing. 

The  encroachment  battleground  involves  debates  over  ill-defined  property  rights, 
a  highly  complex  concept.  Ownership  of  land  is  actually  a  set  of  property  rights 
that  contains  many  components — some  of  which  are  spelled  out  in  legal  docu¬ 
ments  and  some  of  which  are  presumed.  Perhaps  the  most  visible  component  of 
land  ownership  is  a  fenceline  on  a  legal  boundary.  However,  some  property 
rights  may  not  extend  to  the  fenceline  while  others  can  extend  beyond  it.  Noise, 
dust,  and  smoke  impacts  from  training  can  and  do  cross  the  fenceline  legally  be¬ 
cause  a  landowner  has  the  right  to  do  so.  For  example,  smoke  from  a  resident’s 
furnace  or  fireplace  can  leave  the  resident’s  property,  as  can  noise  or  light.  How¬ 
ever,  there  are  limits  set  in  local,  county,  state,  and  Federal  laws  that  must  not 
be  breached. 
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Objectives 

The  objective  of  this  work  was  to  provide  installation  planners  with  better  tools 
to:  (1)  predict  urban  growth  patterns  that  might  occur  in  response  to  proposed 
regional  plans  involving  certain  investments  and  policies,  and  (2)  predict  the 
change  in  training  and  testing  opportunity  that  might  occur  as  a  result  of  those 
patterns. 


Approach 

Prediction  of  urban  patterns  is  accomplished  using  the  Land  Use  Change  (LUC) 
model,  which  is  part  of  a  suite  of  tools  called  the  military  Land  use  Evolution 
and  Assessment  Model  (mLEAM).  Researchers  outlined  the  specific  steps  in¬ 
volved  in  preparing  input  maps  for  the  mLEAM  LUC  model,  with  accompanying 
specific  (case  study)  examples,  and  all  data  and  supporting  programming  scripts. 


Scope 

This  document  provides  information  that  will  assist  GIS  technicians  in  the 
preparation  of  input  maps  for  the  mLEAM  LUC  model.  GIS  technicians  will  find 
general  descriptions  of  the  GIS  processing  steps  that  involve  data  acquisition  to 
final  development  of  mLEAM  input  maps  and  the  description  of  the  specific 
steps  used  to  create  maps  for  the  multi-county  area  surrounding  Fort  Bragg,  NC. 
All  data  and  supporting  scripts  are  provided  with  an  accompanying  CD.  Appen¬ 
dix  A  contains  a  listing  of  the  CD  contents. 


Mode  of  Technology  Transfer 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW)  at 
URL: 

http://www.cecer.army.mil 


ERDC/CERL  TR-04-24 


3 


2  Overview 


Projecting  Urban  Growth  to  Predict  Military  Encroachment 

Assume  that  you  are  a  GIS  technician  assigned  to  predict  urban  patterns  around 
a  military  installation  over  20  to  50  years  for  the  purpose  of  predicting  future 
constraints  on  training  and  testing.  The  mLEAM  model  provides  a  set  of  tools 
and  procedures  for  doing  this  in  an  optimal  manner  using  available  geographic 
information  system  (GIS)  software  and  data.  At  the  core  of  mLEAM  is  an  urban 
growth  model  developed  to  efficiently  and  effectively  predict  urban  footprints  in 
a  manner  suitable  for  predicting  future  constraints  to  training  and  testing. 
There  are  many  urban  growth  models;  each  is  designed  to  answer  a  specific  set 
of  questions  about  that  growth.  It  is  useful  to  understand  that  the  land  use 
change  (LUC)  model  at  the  core  of  mLEAM  does  not  predict  settlement  patterns 
of  different  socioeconomic  groups,  predict  the  development  of  new  roads  and 
highways,  or  accurately  predict  the  specific  location  of  new  growth.  It  does  not 
take  into  account  any  information  about  the  behavior  of  individual  landowners, 
bankers,  politicians,  or  voters. 

The  inputs  to  mLEAM  are  minimal.  This  way,  the  model  provides  general  use¬ 
fulness  in  a  cost-effective  way  to  a  broad  audience.  The  inputs  are: 

•  raster  GIS  maps: 

-  current  land  use 

-  current  roads  and  highways 

-  no  growth  areas 

-  employment  centers  (e.g.,  city  centers)  -  location  and  size 

•  population 

-  Current 

-  Projections 

•  regional  investments  and  policies 

-  Changes  in  road  network 

-  New  highways 

-  New  ramps  to  limited  access  highways 

-  Change  in  no-growth  areas. 


Optional  inputs  involve  a  number  of  model  calibration  parameters. 
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The  mLEAM  LUC  output  is  a  time  series  showing  changes  in  urban  patterns. 
The  LUC  model  predicts  conversion  of  available  undeveloped  land  to  residential, 
commercial,  and  urban  open  space.  Predicted  urban  patterns  are  then  available 
for  processing  to  identify  changes  in  training  and  testing  opportunities  (which 
will  be  described  in  follow-up  technical  reports). 


Figure  1.  Conceptual  diagram  of  complete  mLEAM  process. 


The  complete  conceptual  mLEAM  approach  (Ligure  1)  is  comprised  of  seven 
steps: 

1.  Collect  necessary  data  -  maps  and  population  projections. 

2.  Capture  proposed  investments  and  policies. 

3.  Prepare  data  for  analysis. 

4.  Predict  land  development  probability. 

5.  Identify  constraint  pressures  on  military  training  and  testing. 

6.  Identify  the  value  of  land  for  training  and  testing  (or  carrying  capacity). 

7.  Compare  alternative  regional  policy/investment  proposals  (COAs). 

Steps  1-4  are  done  mostly  on  a  computer  within  an  iterative  framework  de¬ 
scribed  in  Chapter  3.  The  result  of  step  4  will  be  one  or  more  mLEAM  graphical 
outputs.  In  steps  5-7,  planners  analyze  and  interpret  the  results  of  the  mLEAM 
output(s)  from  step  4  to  develop  future  plans  for  land  in  and  around  the  military 
installation.  Each  of  these  seven  steps  is  detailed  in  Chapter  3. 

While  conceptually  simple,  the  procedures  described  in  this  document  for  con¬ 
ducting  land  use  evolution  studies  are  best  implemented  by  GIS  technicians  and 
computer  programmers. 
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The  mLEAM  LUC  Algorithm 

As  discussed  above,  the  LUC  model  urban  growth  projection  takes  a  rather  small 
set  of  user  inputs  and  turns  maps  showing  future  urban  growth  patterns.  The 
technical  approach  taken  to  process  the  inputs  into  the  outputs  is  presented  suc¬ 
cinctly  here  and  is  intended  primarily  for  managers  and  supervisors.  The  three 
gross  steps  are: 

•  acquire  and  prepare  base  maps  (about  2  days) 

•  process  maps  using  GIS  into  inputs  for  the  LUC  model  (about  1  day) 

•  run  the  LUC  model  and  visualize  results  (about  1  day  per  scenario). 

The  first  step  primarily  involves  acquiring  nationally  available  maps  from  any  of 
a  number  of  specific  government  sources.  mLEAM  was  designed  to  be  generally 
useful  and  requiring  only  generally  available  information  meets  this  goal,  al¬ 
though  locally  available  information  can  be  used  to  provide  more  accurate  in¬ 
puts.  Acquisition  and  preparation  of  the  required  maps  can  take  less  than  a  day; 
the  sources  and  maps  are  presented  in  detail  in  later  sections.  This  step  also  re¬ 
quires  hands-on  GIS  processing  to  turn  the  acquired  information  into  specifically 
formulated  maps  so  that  the  automatic  processing  in  the  third  step  can  be  per¬ 
formed.  Required  processing  involves  putting  all  of  the  maps  into  a  common 
equal-area  projection,  giving  the  maps  specific  names,  and  ensuring  that  the 
contents  of  the  maps  use  specific  category  values.  Acquisition  and  processing 
should  take  less  than  2  solid  working  days  with  a  modern  GIS. 

The  second  step  involves  running  scripts  that  invoke  specific  GIS  processing. 
Scripts  are  available  for  users  of  the  ESRI  Arc  Info  software  on  Windows-based 
PCs  and  for  users  of  the  GRASS  GIS  software  on  Unix/Linux  machines.  In  this 
step,  travel  times  are  computed  for  each  cell  (location)  to  such  things  as  the  near¬ 
est  small,  medium,  and  large  cities;  water  and  forest;  and  roads,  highways,  and 
interstates.  This  step  can  take  very  little  operator  time  although  the  computer 
processing  of  a  multi-county  area  can  involve  several  hours.  Inspection  of  results 
is  required  and  may  result  in  the  modification  of  input  maps  and  parameters. 

The  third  step  involves  running  the  LUC  (land  use  change)  model.  In  this  step, 
the  model,  a  configuration  file,  and  the  data  prepared  in  the  third  step  must  be 
brought  together.  The  model  can  take  several  hours  to  run  a  large  area  on  a  sin¬ 
gle  processor  computer.  The  model  is  written  in  C  and  is  generally  run  on  a 
Unix-based  (e.g.,  Linux)  computer.  Output  from  the  model  can  be  taken  back 
into  a  GIS  for  analysis  and  visualization. 
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The  results  of  the  processing  can  then  be  processed  to  understand  the  impacts  of 
the  predicted  patterns  on  the  opportunities  for  training  and  testing.  At  the  time 
of  this  writing,  procedures  were  being  developed  to  automate  this  analysis  and 
will  appear  in  future  reports. 

The  remainder  of  the  information  in  this  document  is  intended  for  GIS  analysts. 
Instructions  identify  how  to  acquire  data  and  software  and  how  to  process  read¬ 
ily  available  data  into  urban  growth  projections  for  your  area  of  interest. 
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3  mLEAM  Procedures 


While  the  conceptual  steps  given  in  Chapter  2  provide  a  theoretical  framework 
by  which  to  envision  the  mLEAM  process,  the  actual  GIS  simulation  deviates 
from  a  strict  step-by-step  method.  In  reality,  the  GIS  (computer)  portion  of  a 
mLEAM  model  is  actually  more  of  an  iterative  process  than  a  step-by-step  proc¬ 
ess,  and  this  is  consistent  with  the  use  of  GIS  as  a  practical  tool.  The  computer 
portion  of  mLEAM  can  be  seen  as  a  separate  but  intertwined  framework  that 
runs  on  top  of  the  conceptual  framework  given  in  Chapter  2. 

The  following  procedures  for  generating  urban  growth  predictions  are  to  be  fol¬ 
lowed  by  individuals  with  significant  experience  in  the  operation  of  geographic 
information  systems.  Today,  GIS  software  makes  different  GIS  functionality 
available  to  everyone  from  novices  to  experts.  These  procedures  involve  the  abil¬ 
ity  to: 

•  Download  maps  from  various  web  sites. 

•  Convert  these  to  maps  with  specific  projections,  names,  and  categories. 

•  Read,  write,  and  execute  involved  scripts. 

Organizing  and  running  the  computer  portion  of  a  mLEAM  model  consists  of  a 
number  of  steps.  While  each  of  these  steps  is  not  inherently  difficult,  they  do  as¬ 
sume  the  user  is  familiar  with  a  number  of  computer  programs.  Instructions  are 
provided  that  assist  in  the  use  of  ESRI’s  suite  of  software  and  the  Unix/Linux 
based  GRASS  software.  Users  can  choose  to  use  one  system  or  the  other,  but  ac¬ 
cess  to  the  strengths  of  both  provides  more  flexibility.  Detailed  instructions  for 
running  the  computer  portion  of  mLEAM  are  included  in  the  Appendices  to  this 
report. 

The  mLEAM  computer  framework,  depicted  in  Figure  2,  consists  of  manipulat¬ 
ing  raw  data  inputs  into  an  overall  land  use  projection,  via  two  main  transforma¬ 
tion  stages.  The  “manipulation”  basically  entails  selecting  germane  information 
from  a  file  and  saving  it  in  a  readable  computer  format.  First,  raw  data  files  are 
collected  from  various  sources  such  as  the  Internet.  Next,  a  GIS  program  such 
as  Arc/Info  or  GRASS  is  used  to  assemble  the  raw  data  files  into  useable  base 
files;  this  is  done  by  hand.  Then  an  automated  GIS  program  takes  the  base  files 
and  manipulates  them  into  Land  Use  Change  (LUC)  input  files. 
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Figure  2.  mLEAM  computer  framework. 


Finally,  the  automated  LEAM  scripts  convert  the  LUC  input  files  into  an  overall 
land  use  projection.  Due  to  the  iterative  nature  of  the  GIS  process,  some  of  the 
intermediate  or  final  maps  may  be  reused  as  starter  or  intermediate  maps  for 
subsequent  LUC  input  files. 

Again,  the  computer  steps  A  through  D  shown  above  are  a  way  to  envision  the 
iterative  nature  of  the  GIS  portion  of  mLEAM,  and  these  steps  are  the  realistic 
way  of  performing  steps  1  through  4  of  the  conceptual  framework  outlined  in 
Chapter  2.  Lor  example,  a  user  may  download  a  raw  data  file  from  the  Web, 
manipulate  it  into  a  base  file  and  then  a  LUC  input  file  to  get  it  ready  for  an  ur¬ 
ban  growth  projection.  The  user  may  then  take  this  file  and  use  it  as  one  of  the 
inputs  for  a  different  base  file.  There  are  more  than  a  dozen  related  LUC  input 
files  to  create,  so  there  are  many  instances  of  re-using  a  common  file  to  create 
new  base  files. 

Still,  the  7-step  process  holds  for  showing  the  mLEAM  procedure  in  a  conceptual, 
sequential  manner.  Ligure  3  shows  these  steps. 
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Figure  3.  Complete  mLEAM  analysis  steps. 


1.  Collect  Necessary  Data 

The  user  will  gather  maps  and  other  data  from  various  sources  (usually  the 
Internet)  for  input  into  a  GIS. 


2.  Capture  Proposed  Investments  and  Policies 

The  user  will  gather  proposed  investments  and  policies  for  the  area,  such  as  new 
highways  or  new  zoning  ordinances — usually  available  only  through  contacts  in 
local  or  state  governing  boards.  Since  these  things  often  affect  the  landscape, 
the  user  will  create  additional  base  maps  that  incorporate  the  proposed  invest¬ 
ments  and  policies.  These  additional  maps  will  act  as  the  basis  for  secondary, 

tertiary,  etc.  scenarios,  so  as  to  answer  questions  like  “If  this  _  is 

built,  how  will  the  new  development  scenario  compare  to  the  ‘business  as  usual’ 
scenario?” 
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Figure  4.  Converting  raw  maps  to  mLEAM  input  maps. 
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3.  Prepare  Data  for  Analysis 

As  Figure  4  shows,  the  user  will  select  pertinent  information  from  the  raw  data 
and  put  it  into  a  computer  format  useable  by  the  mLEAM  software. 

There  are  two  methods  for  creating  the  mLEAM  input  maps:  (1)  using  ESRI  Ar- 
cGIS  exclusively,  or  (2)  using  ESRI  ArcGIS  and  GRASS  GIS.  The  advantage  of 
using  ESRI  ArcGIS  exclusively  is  that  the  user  can  perform  all  operations  with 
one  program  on  Windows.  The  advantage  of  using  ESRI  ArcGIS  in  conjunction 
with  GRASS  is  that  GRASS  has  more  capabilities  for  handling  raster  data.  This 
method  requires  the  use  of  ESRI  ArcGIS  very  briefly  in  the  initial  stages  of  the 
mapmaking  process,  since  the  downloaded  files  usually  come  in  ESRI  Shapefile 
format.  This  method  also  requires  access  to  both  a  Windows  machine  and  a 
UNIX  machine. 

Create  Land  Use  Change  Input  Maps 

The  goal  is  to  create  several  specific  Land  Use  Change  input  maps  (Table  1), 
which  are  described  in  the  following  sections. 


Table  1.  Land  use  change  input  maps. 


Map 

Description 

Boundary 

Boundary  of  study  region  (usually  a  group  of  counties) 

Land  Cover 

Land  cover  of  study  region  (NLCD  data) 

Metrobuffer 

Buffer  around  municipal  boundaries  in  study  region  (typically  1.5  miles) 

No  Growth 

Areas  in  study  region  in  which  development  is  not  possible  or  highly  unlikely 

Slope 

Slope  of  terrain  in  study  region 

Forest  Attractor 

Distance  to  nearest  forest 

Water  Attractor 

Distance  to  nearest  water  body 

Highway  Attractor 

Travel  time  to  nearest  highway 

Intersection  Attractor 

Travel  time  to  nearest  major  road  intersection 

Ramp  Attractor 

Travel  time  to  nearest  highway  ramp 

Road  Attractor 

Travel  time  to  nearest  road 

Cities  Attractor  (large, 

Travel  time  to  nearest  city  (divided  into  three  groups — small,  medium,  and 

medium,  small) 

large) 

Boundary.  This  map  simply  indicates  the  boundaries  of  LEAM  simulation 
runs.  Each  grid  cell  in  which  the  model  will  be  run  is  indicated  with  the  number 
1  and  all  other  cells  are  given  the  value  of  zero. 

Land  Cover.  The  National  Land  Cover  Data  (NLCD)  classification  system  is 
used.  Figure  5  shows  the  raster  map  showing  color-coded  categories  below  use 
the  NLCD  category  definitions. 
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Figure  5.  Sample  land  cover  map. 


The  NLCD  category  numbers  represent  land  cover  as  follows: 

11 .  Open  Water  41 .  Deciduous  Forest 

12.  Perennial  Ice/Snow  42.  Evergreen  Forest 

21.  Low  Intensity  Residential  43.  Mixed  Forest 

22.  High  Intensity  Residential  51.  Shrub  land 

23.  Commercial/Industrial/Transportation  61.  Orchards/Vineyards/Other 

31.  Bare  Rock/Sand/Clay  71.  Grasslands/Herbaceous 

32.  Quarries/Strip  Mines/Gravel  Pits  81.  Pasture/Hay 

33.  Transitional  82.  Row  Crops 

83.  Small  Grains 

84.  Fallow 

85.  Urban/Recreational  Grasses 

Further  details  on  these  definitions  can  be  found  at  URL: 
http://landcover.usqs.gov/classes.html 


Metrobuffer.  Cells  defined  as  within  a  growth  sphere  of  influence  are  indi¬ 
cated  with  a  1,  others  are  zero  (Figure  6).  This  map  can  capture  the  extent  of  the 
anticipated  legal  boundaries  of  cities  over  the  course  of  a  LEAM  simulation. 
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Figure  6.  Sample  metro  buffer  map. 


Nogrowth.  The  nogrowth  map  is  also  a  one-zero  map.  Cells  that  cannot  grow 
are  indicated  with  the  value  of  1.  In  the  sample  map  below,  the  cells  coded  with 
a  “1”  are  associated  with  interstate  highway  rights  of  way.  Other  examples  of 
no-growth  cells  are  those  associated  with  water,  no-growth  zones,  parks,  forests, 
nature  preserves,  etc.  This  map  (Figure  7)  is  used  to  capture  many  alternative 
land  management  policies  (e.g.,  zoning)  and  investments  (e.g.,  purchase  of  prop¬ 
erty  and/or  property  rights.) 
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Figure  7.  Sample  no  growth  map. 
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Slope.  The  slope  map  (Figure  8)  indicates  the  slope  of  the  cell  in  degrees. 

The  slope  map  is  used  in  the  LEAM  process  because  the  likelihood  of  develop¬ 
ment  is  generally  higher  in  areas  with  gentle  slopes  and  generally  lower  in  areas 
with  steep  slopes. 
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Forest_attractor.  The  forest  attractor  is  simply  the  straight-line  distance  to 
the  nearest  cell  identified  as  having  forest  in  the  land  cover  map.  Cells  contain¬ 
ing  forest  are  coded  as  “0”  (zero  distance  to  the  nearest  forest).  The  cell  size  the 
map  shown  in  Figure  9  is  30-meters  and  the  distances  to  nearby  cells  containing 
forest  is  obviously  recognizing  this  cell  size.  The  forest  attractor  map  is  used  in 
LEAM  because  the  existence  of  trees  makes  development  more  attractive. 
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Figure  9.  Sample  forest  attractor  map. 


Water_attractor.  The  water  attractor  map  is  identical  in  concept  to  the  forest 
attractor.  The  sample  map  portion  shown  in  Figure  10  indicates  a  body  of  water 
above  the  area  just  to  the  left  of  center.  The  water  attractor  map  is  used  in 
LEAM  because  the  existence  of  water  makes  development  more  attractive. 
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Figure  10.  Sample  water  attractor  map. 
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Highway_attractor.  This  and  the  following  three  maps  identify  the  driving 
time  in  minutes  to  the  nearest  state  highway,  intersections  of  open  roads,  lim¬ 
ited-access  highway  ramps,  and  main  roads.  The  driving  time  is  calculated  with 
an  algorithm  that  takes  into  account  speed  limits  for  roads,  as  defined  by  the 
user.  The  creation  of  these  maps  relies  on  the  development  of  a  cell-traversal 
travel  time  map.  Then,  GIS  analyses  that  identify  the  minimum  travel  time  to 
the  respective  resources  are  conducted  to  create  these  four  maps.  In  the  images 
below,  black  represents  zero  travel  time  and  the  whiter  the  cell,  the  greater  the 
travel  time.  The  pure  white  areas  indicate  no-growth  locations  (water,  Federal 
lands,  etc.)  upon  which  travel  is  near  impossible;  hence,  travel  times  are  arbi¬ 
trarily  set  extraordinarily  high.  The  highway  attractor  map  (Figure  11)  is  used 
in  LEAM  because  development  is  more  likely  near  a  highway. 


Figure  11.  Sample  highway  attractor  map. 


lntersection_attractor.  The  map  shown  in  Figure  12  depicts  the  travel  time  to 
the  nearest  highway  intersection.  The  intersection  attractor  map  is  used  in 
LEAM  because  development  is  more  likely  near  highway  intersections. 


Figure  12.  Sample  intersection  attractor  map. 
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Ramp_attractor.  The  map  shown  in  Figure  13  depicts  the  travel  time  to  the 
nearest  limited  access  highway  ramp  (e.g.,  U.S.  Interstates).  The  intersection 
attractor  map  is  used  in  LEAM  because  development  is  more  likely  near  ramps. 


Figure  13.  Sample  limited  access  highway  ramp  attractor  map. 

Road  attractor.  The  map  shown  in  Figure  14  depicts  the  travel  time  to  the 
nearest  road.  The  road  attractor  map  is  used  in  LEAM  because  the  existence  of 
a  road  increases  the  likelihood  of  development  occurring. 


Figure  14.  Sample  county  road  attractor  map. 

Cities_attractor.  The  cities  attractor  map  (Figure  15)  denotes  driving  time  to 
the  nearest  employment  centers  (municipalities).  Because  there  are  many  mu¬ 
nicipalities  in  a  study  area,  they  are  grouped  into  a  few  categories  by  population 
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(e.g.,  large,  medium  and  small  cities).  Attractor  map  is  created  for  each  group  of 
municipalities  by  the  same  way  as  road  attractor  maps.  Attractor  map  for  each 
group  is  given  different  weights. 


General  Approach 

The  mLEAM  modeler  must  develop  the  input  maps  (shown  above)  for  their  area 
of  interest.  There  are  many  approaches  that  can  be  used  and  users  are  encour¬ 
aged  to  explore  options  and  opportunities.  However,  the  following  figure  sug¬ 
gests  the  maps  that  might  be  used  to  create  the  mLEAM  input  maps  (identified 
along  the  right  side).  The  grey-background  maps  along  the  left  side  are  locally 
specific  user  maps  that  capture  plans  that  mLEAM  will  test  with  respect  to  an¬ 
ticipated  urban  growth.  Those  without  the  grey  background  represent  the  ge¬ 
neric  nationally  available  maps.  Two  temporary  maps  (bottom  center)  are  sug¬ 
gested  as  important  steps  in  the  process  of  creating  attractiveness  maps  that 
consider  driving  time  on  the  current  and  planned  road/highway  network.  One  is 
simply  a  description  of  the  driving  time  required  to  cross  each  30-meter-square 
cell  and  is  used  to  generate  several  needed  images.  The  other  indicates  the  fun- 
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damental  economic  draw  of  cities  and  their  associated  job  and  commercial  cen¬ 
ters. 

Most  of  the  map  processing  is  straightforward.  Creation  of  the  slope  map  re¬ 
quires  running  standard  slope-aspect  analysis  models  available  in  most  raster 
GIS  packages.  The  boundary  map  is  simply  an  overlay  identifying  the  counties 
involved  in  the  simulation  process.  Similarly,  the  Metrobuffer  map  is  a  1-0  map 
identifying  the  extent  of  anticipated  edges  of  cities  through  the  course  of  a  simu¬ 
lation.  These  and  the  no-growth  map  are  most  challenging  with  respect  to  ex¬ 
tracting  agreement  from  the  user  community.  Forest  and  water  attractiveness 
are  simply  proximity  analyses  that  give  distance  from  each  cell  to  the  closest  cell 
containing  forest  or  water. 

The  road,  highway,  ramp,  and  intersection  attractiveness  maps  are  based  on  the 
travel  time  to  drive  to  these  features.  A  travel  time  map  is  first  developed, 
which  identifies  the  time  (in  minutes)  required  to  traverse  each  cell.  Artificially 
high  values  are  assigned  to  such  areas  as  water,  permanent  easements  for  lim¬ 
ited  access  roads/highways,  and  railroads,  to  essentially  indicate  that  passage  is 
impossible.  Using  a  raster  GIS  cumulative  cost  function,  this  map  can  be  used  to 
generate  the  various  travel-time  based  attractiveness  maps. 

The  most  challenging,  and  perhaps  the  most  important,  input  map  to  generate  is 
the  cities_attractor  map.  This  map  captures  the  notion  of  travel  time  to  jobs  and 
shopping.  In  the  current  release  a  value  is  assigned  to  each  cell  with  respect  to 
the  most  attractive  center  (rather  than  a  cumulative  summing  of  the  attractive¬ 
ness  to  all  nearby  centers).  “Attractiveness”  in  this  case  is  defined  as  the  city 
center  closest  to  a  grid  cell  with  respect  to  shortest  travel  time. 


4a.  Predict  Land  Development  Probability,  or  4b.  Predict  a  Number  of 
Alternative  Futures 

The  user  will  run  the  mLEAM  simulation(s)  with  the  mLEAM  input  maps  cre¬ 
ated  in  step  2,  and  one  or  more  mLEAM  output  scenarios  will  result.  The  output 
map  will  show  the  likelihood  of  urban  development  for  a  given  region.  If  more 
than  one  input  scenarios  are  run,  the  user  will  be  able  to  compare  the  outputs  for 
multiple  scenarios,  thus  allowing  the  prediction  of  a  number  of  alternative  de¬ 
velopment  futures. 


ERDC/CERL  TR-04-24 


19 


5.  Identify  Constraint  Pressures  on  Military  Training  and  Testing 

Since  the  mLEAM  outputs  are  purely  graphical,  planners  can  easily  glean  trends 
from  the  mLEAM  models  and  see  which  areas  of  an  installation  might  be  af¬ 
fected  by  urban  encroachment. 


6.  Identify  the  Value  of  Land  for  Training  and  Testing  (or  Carrying 
Capacity) 

Military  planners  must  decide  the  comparative  values  of  installations’  lands  with 
respect  to  training  and  testing.  In  many  cases,  training/testing  operations  on 
certain  areas  within  an  installation  have  been  suspended  due  to  massive  urban 
encroachment  and  the  noise  and  dust  complaints  and  lawsuits  that  follow.  In 
these  cases,  it  has  simply  become  too  costly  for  the  military  to  use  those  grounds 
for  training  and  testing.  It  is  incumbent  upon  the  military  planners  to  recognize 
which  lands  are  of  top  priority  for  saving  and  which  are  relatively  expendable. 


7.  Compare  Alternative  Regional  Policy/Investment  Proposals  (COAs) 

The  final  goal  of  the  mLEAM  approach  is  for  planners  to  take  the  insights 
learned  from  the  growth  models  and  prepare  regional  policies  that  will  benefit 
the  region  the  most.  The  graphical  models  can  help  show  the  predicted  end  re¬ 
sult  of  a  particular  investment  proposal  or  policy,  and  the  opportunity  to  model 
multiple  scenarios  allows  for  easy  comparison  of  alternative  plans.  If  one  of  the 
models  predicts  growth  in  a  pattern  that  is  preferable  to  the  others,  planners 
might  use  that  regional  policy  (or  aspects  of  it)  as  a  starting  point  for  shaping 
urban  growth  in  their  area. 
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4  An  mLEAM  Example:  Fort  Bragg,  NC 


Challenge 

Fort  Bragg,  a  large  Army  base  in  south-central  North  Carolina,  has  been  experi¬ 
encing  urban  encroachment  on  its  borders  from  neighboring  Fayetteville,  result¬ 
ing  in  more  noise  and  dust  complaints  from  the  new  residents.  This  caused  the 
installation  to  reconsider  its  use  of  some  border  lands  for  training  and  testing. 
The  loss  of  this  land  for  training  and  testing  has  both  economic  and  national  se¬ 
curity  implications.  Due  to  the  construction  of  a  new  four-lane  limited-access 
highway  running  between  Fort  Bragg  and  Fayetteville,  urban  encroachment  was 
foreseen  as  only  getting  more  intrusive.  The  challenge  was  for  a  model  to  accu¬ 
rately  and  precisely  predict  where  the  new  growth  would  occur,  both  with  the 
new  highway  and  without  it. 


Collecting  Necessary  Data  and  Preparing  for  Analysis 

This  section  will  describe  the  process  by  which  the  Fort  Bragg  area  data  was  col¬ 
lected  and  converted  into  useable  base  files.  This  corresponds  with  the  steps  A, 
B,  and  C  of  the  mLEAM  computer  approach,  described  in  the  first  section  of 
Chapter  3.  Again,  since  GIS  is  an  iterative  process,  there  will  not  be  a  strict  ad¬ 
herence  to  following  a  sequence  from  A  to  B  to  C.  This  section  will  show  both 
methods  for  producing  the  mLEAM  input  maps:  (1)  using  ESRI  Arc/Info  exclu¬ 
sively,  and  (2)  using  ESRI  ArcGIS  alongside  GRASS. 

Following  the  section  describing  the  steps  taken  in  creating  the  input  maps,  a  set 
of  schematic  diagrams  is  included.  These  schematic  diagrams  depict  the  order  of 
data  acquisition  and  manipulation  for  each  map  layer.  There  are  four  diagrams: 

1.  Scenario  without  new  highway,  done  with  ESRI  ArcGIS  exclusively. 

2.  Scenario  without  new  highway,  done  with  ESRI  and  GRASS. 

3.  Scenario  with  new  highway,  done  with  ESRI  ArcGIS  exclusively. 

4.  Scenario  with  new  highway,  done  with  ESRI  and  GRASS. 

The  term  “BR”  is  generally  used  as  an  abbreviation  for  “Bragg”  (to  shorten  the 
filenames). 
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CD  of  Files  Used  to  Create  Fort  Bragg  mLEAM  Scenario 

A  CD  was  burned  containing  the  files  used  in  the  creation  of  the  Fort  Bragg 
mLEAM  scenario.  In  most  cases,  the  CD  filenames  correspond  with  the  sug¬ 
gested  filenames  given  in  the  tutorial  shown  on  the  following  pages.  A  complete 
list  of  the  files  on  the  CD  is  included  in  the  Appendix  A. 


ESRI  ArcGIS  Only  Version 
Initial  Setup 

1.  Create  directory  on  hard  drive  (e.g.,  C:\Fort_Bragg)  to  hold  data.  There  should 
not  be  any  spaces  in  the  path  names. 

2.  Obtain  Arc/Info  scripts  and  support  files  from  LEAM  lab,  and  save  them  to  the 
directory  created  in  step  1.  These  files  are: 


1bnd_cover 

12slope 

geog_alb 

1  bnd_g  rid 

1 3stateatr 

growremap 

2landuse 

14treeatr 

intersection_point 

3munbuf 

1 5wateratr 

munbndremap 

4travel 

16facility 

nogrowremap 

6roadatr 

astrgrdp 

roadremap 

7flood 

astrgrdp2 

travelremapl 

8cityatr 

boundremap 

travelremap2 

8growth 

boundremap2 

travelremap3 

9nogrow 

cityremap 

travelremap4 

lOrampatr 

11 rdintatr 

facilityremap 
fid  remap 

travelremap5 

3.  Create  workspace  in  Arc/Info.  The  command  will  be  “w”  followed  by  the  path  to 
the  directory  created  in  step  1  (e.g.,  “Arc:  w  c:\fort_Bragg”).  This  workspace  will 
be  the  place  where  files  are  accessed  and  saved.  This  workspace  must  be  opened 
each  time  Arc/Info  is  used. 

4.  Notes:  Some  of  the  following  maps  can  be  created  by  running  AML  scripts.* 

These  scripts  are  collections  of  commands  that  run  automatically  in  Arc/Info.  If 
the  scripts  are  used,  some  of  the  intermediate  steps  described  will  not  be  neces¬ 
sary,  because  the  scripts  already  contain  those  commands.  The  references  to  the 
scripts  are  shown  in  italics. 

It  is  imperative  that  the  AML  scripts  be  saved  to  the  same  directory  as  the  maps 
being  worked  on;  otherwise,  Arc/Info  will  not  recognize  them.  When  running  an 


Appendix  B  includes  a  listing  of  AML  scripts. 
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AML  script  from  an  “Arc:”  prompt,  type  “&run  <pathname>.<AMLscript>  Ar¬ 
guments  (if  necessary)>”  without  the  quotes. 

Creating  Maps 

1.  County  Boundary  Map.  This  will  describe  taking  a  raw  vector  data  file 
of  county  boundaries  from  a  data  source  and  manipulating  it  into  a  LEAM  input 
map  file.  The  result  will  be  two  files:  (1)  a  vector  coverage  of  the  county  bounda¬ 
ries  in  the  Bragg  area,  and  (2)  a  grid  of  the  same  area. 

Source(s) 

http://www.census.qov/qeo/www/cob/co2000.html 

Procedure 

1.  Go  to  website  above,  select  North  Carolina,  download  and  unzip  the  file 
co37_dOO.eOO. 

2.  In  Arc/Info,  Import  the  cover  by  Arc  IMPORT  command  (e.g.,  “Arc:  import  cover 
co37_dOO.eOO  nc_counties”). 

3.  Project  NC_COUNTIES  using  geog_alb.txt  projection  file  (e.g.,  “Arc:  project  cover 
nc_counties  nc_counties_al  geog_alb.txt”). 

4.  Build  polygons  by  CLEAN  command  (e.g.,  “Arc:  clean  nc_counties_al”)- 

5.  To  limit  the  area  to  just  the  Fort  Bragg  area,  select  necessary  counties  using  the 
RESELECT  command  (Counties  in  Fort  Bragg  region:  Cumberland,  Harnett, 
Hoke,  Lee,  Moore,  Richmond,  and  Scotland).  (“Arc:  reselect  nc_counties_al 
br_county  poly”  “  >:  RES  county  lk  ‘051’  or  county  lk  ‘085’  or  county  lk  ‘093’  or 
county  lk  ‘105’  or  county  lk  ‘125’  or  county  lk  ‘153’  or  county  lk  ‘165’”  |  type  ‘n’s  to 
the  questions  that  follow. 

AML  Script  Opportunity:  Procedures  (2)  thru  (5)  can  be  done  by  running  the 
boundary  cover  AML  script  (e.g.,  “Arc:  &run  C:  \Fort_Bragg  l  lbndjcover.aml”): 

-  For  the  exchange  file  name,  type  “CO37JD00”. 

-  For  the  resulting  cover  name,  type  “BR_COUNTY”. 

-  For  the  logical  expressions,  type  the  commands  I  arguments  used  in  step  (5) 
at  the  >:  prompt. 

6.  Convert  BR_COUNTY  to  a  grid  (e.g.,  BR_COUNTY_G)  and  reclassify  so  that  all 
cells  have  the  value  of  1  (e.g.,  BR_GRID_1). 

7.  Clip  BR_GRID_1  with  the  cover  BR_COUNTY  as  BR_GRID  using  GRIDCLIP 
command  in  GRID  module  (e.g.,  “Grid:  gridclip  br_grid_l  br_grid  cover 
br_county”). 

8.  Export  BR_COUNTY  as  an  eOO  file  so  that  it  can  be  used  for  other  applications 
(e.g.,  “Arc:  export  cover  br_county  br_county.eOO”). 
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9.  Export  BR_GRID_G_1  as  ASCII  (e.g.,  “Arc:  gridascii  br_grid_g_l  br_county.asc”). 

AML  Script  Opportunity:  Procedures  (6)  through  (9)  can  be  done  by  running  the 
boundary  grid  AML  script:  “Arc:  &run  C:  \Fort_Bragg  l  lbnd_grid.aml”: 

-  For  boundary  cover  name,  type  “BRjCOUNTY”. 

-  For  resulting  grid  name,  type  “BR_GRID”. 
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2.  Landuse  Map.  This  will  describe  taking  a  raw  NLCD  raster  data  file 
from  a  website  and  manipulating  it  into  a  LEAM  input  map  file.  The  result  will 
be  a  raster  map  showing  land  use  in  the  Bragg  area  by  grid  cell  (Figure  16). 


Figure  16.  Land  use,  Fort  Bragg  region. 


Sourcefs) 

USGS  National  Land  Cover  Data  (NLCD)  from: 

http://edcsqs9.cr.usqs.qov/pub/data/landcover/states/ 

Road  network  data  from  Census: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 

Municipal  boundary  data  from  Census: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 

Procedure 

1.  Go  to  USGS  NLCD  website,  Choose  North  Carolina,  Download  a  NLCD  TIFF  file 
of  North  Carolina  from  NLCD  site,  unzip  and  save  as  grid  NC_NLCD. 

2.  Download  road  shapefiles  and  municipal  boundary  shapefiles  of  the  same  coun¬ 
ties  (Cumberland,  Harnett,  Hoke,  Lee,  Moore,  Richmond,  Scotland)  and  merge 
into  one  file  for  each  (e.g.,  BR_ROAD.SHP  and  BR_MUNBND.SHP). 

3.  To  do  this:  From  ESRI  website:  choose  North  Carolina  ->  Select  by  layer  “Line 

Features — Roads”  Pick  counties  Download  files  and  unzip  Open  layers  in 

ArcMap  and  merge  into  one  road  shapefile 

Repeat  for  municipal  boundaries,  except  choose  “Designated  Places  2000”  from 
the  Select  by  Layer  list. 

4.  In  Arc/Info,  clip  the  NC_NLCD  grid  with  the  boundary  cover  BR_COUNTY 
save  as  BR_NLCD  (e.g.,  “Arc:  grid”  “Grid:  gridclip  nc_nlcd  br_nlcd  cover 
br_county”). 
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5.  Convert  BR_ROAD.SHP  to  a  cover  using  SHAPEARC  and  ROUTEARC  com¬ 
mands  in  ARC  (e.g.,  “Arc:  shapearc  br_road.shp  br_roadl  subclass”  ->  “Arc:  rou- 
tearc  br_roadl  subclass  br_road2”)  then  kill  BR_ROADl  (“Arc:  kill  br_roadl  all”). 

6.  Project  the  cover  BR_ROAD2  to  Albers  as  BR_ROADl.  (This  step  may  be  unnec¬ 
essary  if  the  road  file  is  already  projected.) 

7.  Convert  BRJV1UNBND.SHP  to  a  polygon  cover  using  SHAPEARC  and  CLEAN 
commands  in  ARC  and  project  to  Albers  as  BRJVIUNBNDl  (“Arc:  shapearc 
br_munbnd.shp  br_munbnd”  ->  “Arc:  clean  br_munbnd”  ->  “Arc:  projectdefine 
Project. 

8.  Select  CFCC2  Al,  A2,  A3  and  A6  class  roads  from  BR_ROADl  and  save  as 
BR_ROAD2  (e.g.,  “Arc:  RESELECT  br_roadl  br_road2  LINE”  then  type  logical 
expressions  like  “res  cfcc2  =  Al'  OR  cfcc2  =  'A2'  OR  cfcc2  =  A3'  OR  cfcc2  =  'A6'”) 
then  kill  BR_ROADl.  The  Al’  and  A2’,  etc.  part  is  case  sensitive. 

This  selection  is  done  to  select  the  major  roads  (Interstates,  Highways,  Ramps). 

9.  Convert  BR_ROAD2  to  a  grid  BR_ROADl  and  reclassify  as  24  (e.g.,  “Grid: 
br_road24  =  (br_roadl  *  0)  +  24”)  then  kill  BR_ROADl. 

10.  Create  70m  buffer  along  BR_ROAD2  as  BR_R_BUFF  and  kill  BR_ROAD2. 

11.  Clip  the  buffer  cover  BR_R_BUFF  using  the  municipal  boundary  cover 
BR_MUNBND1  as  BR_R_BUFF1  (e.g.,  “Arc:  ERASE  br_r Jbuff  br  munbnd  1 
br_r_huffl”)  and  clean  BR_R_BUFF1. 

12.  Convert  BR_R_BUFF1  to  a  grid  BR_R_BUFF2  with  INSIDE  field  as  cell  values 
and  kill  BR_R_BUFF1. 

13.  Select  cell  value  100  from  BR_R_BUFF2  and  save  as  another  grid  BR_R_BUFF1. 

14.  Select  cell  values  21  and  22  from  BR_NLCD 

(21  =  Low  Intensity  Residential,  22  =  High  Intensity  Residential). 

AML  Script  Opportunity:  Procedures  (4)  through  (13)  can  be  done  using  the  lan- 
duseAML  script:  “Arc:  &r  C:\Fort_Bragg\2landuse_new.aml  br_road  br_nlcd 
brjmunbnd  br_nlcd_r” 

Caution:  There  are  four  arguments  used  in  this  AML  script:  The  first 
argument  is  the  name  of  road  shape  file,  the  second  one  is  the  name  of 
land  cover  grid,  the  third  one  is  the  name  of  municipal  boundary 
shapefile  and  the  last  one  is  the  output  grid  name  A  no  extensions 
(SHP). 

Shapefiles  should  be  in  Geographic  Coordinates. 

Remap  files  ASTRGRDRP.  TXT  and  ASTRGRDRP2.  TXT  must  be  lo¬ 
cated  in  the  workspace. 
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3.  Buffer  Map.  This  will  describe  taking  a  raw  vector  data  file  of  munici¬ 
pal  boundaries  from  a  data  source,  combining  it  with  a  county  boundary  file,  and 
manipulating  into  a  LEAM  input  map  file.  The  result  will  be  a  raster  map  show¬ 
ing  a  1.5  mile  buffer  around  all  municipal  boundaries  in  the  Bragg  area  (Figure 
17). 


Figure  17.  Fort  Bragg  municipal  boundary  buffer  (1.5  Mi). 

Sources 

Municipal  boundary  data  from  http://arcdata.esri.com/data/tiqer2000/tiqer  download.cfm 
County  boundary  data  from  #1  (e.g.,  br_county.eOO) 

Procedure 

1.  Go  to  the  web  page  above  and  select  “North  Carolina.” 

2.  Download  designated  places  2000  for  counties  of  interest  (Cumberland,  Harnett, 
Hoke,  Lee,  Moore,  Richmond,  Scotland). 

3.  Unzip  the  file,  open  in  ArcMap  and  merge  into  one  shape  file  (e.g., 
BR_MUNBND.SHP). 

4.  Import  the  boundary  cover  from  the  boundary  data  (e.g.,  “Arc:  import  cover 
br_county.eOO  br_bnd”)  and  convert  to  a  grid  (e.g.,  “Arc:  polygrid  br_bnd 
br_bnd_g”  and  type  ‘30’  for  cell  size  and  ‘y’  for  converting  the  entire  cover). 

5.  Convert  the  shapefile  to  a  coverage  (e.g.,  “Arc:  shapearc  br_munbnd.shp  br_mun 
subclass”  “Arc:  clean  br_mun  br_mun_poly”). 

Note:  The  name  SUBCLASS  is  the  name  of  info  file  where  the  attributes 
are  saved.  It  does  not  really  matter  how  you  name  it. 

6.  Project  the  BR_MUN_POLY  from  GCS  to  Albers,  name  as  BRJV1UN,  and  build 
(or  clean)  BR__MUN. 

Note:  Depending  on  the  input  file  coordinate  system,  it  may  be  necessary 
to  use  PROJECTDEFINE  before  projecting. 

7.  Create  1.5  mile  buffers  outside  the  boundary  of  BR_MUN. 

8.  Create  a  new  field  in  the  table  and  put  1  to  all  the  records  (as  integer)  Arc- 
View-^  Attribute  Tabled  Editing-^  Add  Field. 
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9.  Convert  to  a  grid  with  the  new  field  as  values  BR_MUNBUFF. 

10.  Merge  BRJMUNBUFF  onto  BR_BND_G  and  save  as  BR_MB. 

11.  Kill  BR_MUNBUFF. 

12.  Clip  the  grid  with  the  boundary  cover  imported  from  br_county.e00  (e.g.,  “Grid: 
latticeclip  br_mb  br_bnd  br_munbuff ’)  ->  BRJMUNBUFF  grid  is  the  final  result. 

AML  Script  Opportunity:  Procedures  (4)  thru  (12)  can  be  done  using  the  munici¬ 
pal  buffer  AML  script:  “Arc:  &r  C:  \Fort_Bragg  1 3munbuf.aml. 

Files  ‘boundremap2’  file  and  ‘munbndremap’  must  be  in  the  workspace. 

For  the  municipal  boundary  shapefile  without  extension,  type 
“BR_MUNBND 

For  the  boundary  GRID,  type  “BRjGRID”. 

For  the  BOUNDARY  cover,  type  “BR_CO UNTY”. 
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4.  Travel  Minute  Map.  This  will  describe  combining  the  road  network  file 
with  the  land  use  file  from  previous  steps.  The  result  will  be  a  raster  file  show¬ 
ing  travel  time  across  each  grid  cell  in  minutes,  for  the  Bragg  area  (Figure  18). 
There  is  a  maximum  travel  time  of  approximately  37  minutes  across  a  grid  cell 
due  to  a  programmer-defined  mathematical  limitation.  This  map  will  be  an  in¬ 
termediate  map — it  will  be  used  in  conjunction  with  other  landscape  maps  to 
make  attractor  (travel  time)  maps,  but  will  not  be  used  as  a  land  use  change  in¬ 
put  map  per  se. 


Figure  18.  Travel  time  across  grid  cells,  Fort  Bragg  region. 

Source 

BR_NLCD  map:  Road  network  data  from: 

http://arcdata.esri.com/data/tiger2Q00/tiger  download. cfm 

Procedure 

1.  Download  road  data  from  the  website  above  for  North  Carolina  counties  of  inter¬ 
est  (Cumberland,  Harnett,  Hoke,  Lee,  Moore,  Richmond,  Scotland). 

2.  Open  these  files  in  ArcMap  and  merge  the  shapefiles  into  one  file  (e.g., 
BR_ROAD.SHP).  It  is  acceptable  to  use  the  same  file  from  the  LANDUSE  map 
#2. 

3.  Select  roads  whose  CFCC2  is  Al,  A2,  A3,  A4,  A5,  A6  or  A7. 

4.  Add  a  new  numeric  field  named  ‘class’  and  assign  1  for  Al  records,  2  for  A2  re¬ 
cords,  3  for  A3  records,  etc. 

Editor->Add  field  called  “Newfield”-> Calculate  as  [Cfcc2].Right(l). 

Then,  Add  another  field  called  “class,”  making  sure  that  it  is  categorized 
as  a  number^  Calculate  as  [Newfield].AsNumber. 

Then  delete  the  field  “newfield”. 

Export  this  to  a  .eOO  file. 

5.  Assigning  travel  speed  to  road  data,  Add  a  new  field  and  assign  values  (travel 
speed  in  O.OOlmph)  as  follows: 
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Al:  70000 
A2:  50000 
A3:  40000 
A4:  30000 
A5:  20000 
A6:  30000 
A7:  10000. 

6.  Convert  to  a  grid  (e.g.,  BR_ROAD)  with  the  values  above  as  cell  values. 

7.  Travel  speed  grid. 

Assign  travel  speed  (O.OOlmph)  values  to  different  land  cover  cells  (name 
it  temporarily  XXX). 

NLCD  codes  21,  22  or  23:  10000  (codes  represent  urban  areas:  21=Low 
Intensity  Residential,  22=High  Intensity  Residential,  23=Comm/Ind 
Transportation) . 

Non-urban  land  cover:  500 

Limited-access  hwy  buffer  cells:  1 

Open  water  and  wetlands  cells:  1 

8.  Merge  BR_ROAD  onto  the  XXX  grid  and  save  as  BR_TRAVELSPEED. 

9.  Travel  minute  grid. 

A  travel  minute  grid  is  made  using  the  equation  under  ArcGrid  (e.g., 

“Grid:  TRAVELMIN  =  1000  /  (26.8224  *  BR_SPEEDLIMIT)”). 

Example:  when  speed  limit  is  60mph,  travel  minute  =  1000  (26.8224  * 

60000)  =  0.000621min/m  =  0.621min/km  =  lmin/mile. 

26.8224  is  a  conversion  factor. 

AML  Script  Opportunity:  Procedures  (5)  thru  (9)  can  be  done  using  the  travel 
map  AML  script,  “Arc:  &r  C:  \Fort_Bragg  1 4travel.aml”. 

Note:  the  ROAD  COVER  file  should  he  a  .e00  file.  R  should  have  the 
“Class”  field  from  step  (1).  It  should  be  in  Geog  Coord  System  (like  the 
original  shapefile  was  in),  not  in  a  Projected  System  like  Albers.  If  the  file 
is  in  a  projected  system,  it  will  result  in  an  error. 

For  the  value  item  field  name,  type  “CLASS”. 

For  limited-access  hwy  class,  type  “1”. 

For  ramp  class,  type  “6”. 

Files  ‘travelremap  1’,  ‘travelremap2’,  ‘travelremap3’,  ‘travelremap4’ ,  and 
‘travelremap  5’  should  be  in  the  same  workspace 

10.  Convert  BR_ROAD.SHP  as  a  cover  and  export  as  an  eOO  file  for  future  use,  e.g.,: 

“Arc:  shapearc  br_road.shp  br_road  subclass” 

“Arc:  routearc  br_road  subclass  br_roadl” 

“Arc:  export  cover  br_roadl  br_road.e00” 
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5.  City  Attractor  Maps  (a.  Large,  b.  Medium,  c.  Small).  This  will  describe 
taking  a  raw  vector  data  file  of  municipal  boundaries  in  the  Bragg  area  from  a 
data  source,  adding  population  characteristics  from  the  U.S.  Census  Bureau,  se¬ 
lecting  certain  polygons  based  on  population,  doing  a  cost  calculation  based  on 
the  travel  time  map  made  previously,  and  manipulating  it  into  three  Land  Use 
Change  input  map  files.  The  result  will  be  three  files:  raster  files  showing  travel 
time  from  each  grid  cell  to:  (1)  large  cities,  (2)  medium  cities,  and  (3)  small  cities 
(Figure  19). 


Figure  19.  Travel  time  to  large  city,  Fort  Bragg  region. 

Source 

Municipal  boundary  (Census  population  data  of  incorporated  places)  from: 
http://arcdata.esri.com/data/tiqer2Q00/tiqer  download. cfm 

Procedure 

1.  Download  municipal  boundary  data  or  use  file  from  #3. 

2.  Open  the  table  of  the  shapefile,  add  a  column  “population”  and  input  the  values 
from  Census. 

3.  Divide  the  cities  into  necessary  groups  (e.g.,  large,  medium  and  small)  based  on 
population  distribution.  This  is  a  user-defined  arbitrary  classification. 

4.  Generate  shapefiles  of  large  cities,  medium  cities  and  small  cities  (e.g., 
BRLCITY.SHP,  BR_MCITY.SHP  and  BR_SCITY.SHP). 

5.  Convert  each  shapefile  to  grids  (e.g.,  BR_LCITY,  etc.). 

6.  Run  COSTDISTANCE  command  with  BR_LCITY  as  source  and  TRAVELMIN 
as  cost  with  the  maximum  value  120  (e.g.,  “Grid:  br_lcity_att  =  int  (costdis- 
tance(br_lcity,travelmin))”) . 

7.  Merge  BR_LCITY_ATT  onto  a  boundary  grid  whose  cell  values  are  all  120. 

8.  Repeat  for  medium  and  small  cities  city  attractor  maps  are  produced  for  each 
group  of  cities. 
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AML  Script  Opportunity:  Procedures  (5)  thru  (7)  can  be  done  using  the  city  at¬ 
tractor  script:  “Arc:  &r  C:  l  Fort_Bragg  1 8cityatr.aml” 

For  the  city  shape  file  name,  type  the  name  of  the  city  file  (e.g., 

BR_LCITY.  SHP) 

Repeat  for  other  size  cities  (there  should  be  a  Large,  Medium,  and  Small 
size  city  attractiveness  file  when  complete) 
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6.  Road  Attractor  Map.  This  will  describe  taking  a  raw  vector  data  file  of 
roads  in  the  Bragg  area  from  a  data  source,  doing  a  cost  calculation  based  on  the 
travel  time  map  made  previously,  and  manipulating  it  into  a  LEAM  input  map 
file.  The  result  will  be  a  raster  file  showing  travel  time  from  each  grid  cell  to  the 
nearest  city  road  (Figure  20). 


Figure  20.  Travel  time  to  nearest  road,  Fort  Bragg  region. 

Source 

Road  network  data  from: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 
Travel  minute  grid  from  4  (e.g.,  TRAVELMIN ) 

Procedure 

1.  Import  the  exchange  file  (eOO)  from  4  to  a  cover  (e.g.,  BR_ROAD). 

2.  Select  CFCC2  =  “A4”  AND  Fetype  =  “RD”  and  convert  to  a  grid  (e.g.,  BR_A4RD). 

Note:  “A4”  is  selected  to  highlight  the  city  roads  (class  A4)  only 

3.  Project  the  grid. 

4.  Run  COSTDISTANCE  command  with  BR_A4RD  as  source  and  TRAVELMIN  as 
cost  and  set  the  values  to  integer. 

AML  Script  Opportunity:  Procedures  (1)  thru  (4)  can  be  done  using  the  road  at¬ 
tractor  AML  script:  “Arc:  &r  C:  \Fort_Bragg  1 6roadatr .ami” 

For  the  value  item  field  name,  type  “CLASS”. 

For  roads  and  streets,  type  “4”. 
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7.  No  Growth  Zone  Map.  This  will  describe  taking  raw  vector  data  files  of 
no  growth  areas  in  the  Bragg  area  from  a  data  source,  and  manipulating  it  into  a 
LEAM  input  map  file.  These  areas  consist  of  Federal  lands,  Indian  lands,  water 
bodies,  roads,  landmark  polygons,  etc.  The  result  will  be  a  raster  file  showing  all 
places  within  the  Bragg  area  on  which  development  cannot  occur  (Figure  21). 


Figure  21.  No  growth  areas,  Fort  Bragg  region. 


Source 

1.  National  Atlas  for  Federal  and  Indian  land: 
http://www.nationalatlas.gov 

2.  Other  local  data  about  ownership  (if  possible) 

3.  Landmark  polygon  data  from: 

http://arcdata.esri.com/data/tiqer200Q/tiqer  download. cfm 

4.  Road  network  data  exchange  file  (eOO) 

5.  Boundary  cover  (e.g.,  BR_COUNTY)  and  boundary  grid  (e.g.,  BR_GRID) 

Procedure 

1.  Merge  all  shapefiles  except  roads  and  boundary  into  one  shapefile  (e.g., 
BR_NOGROW). 

2.  Convert  the  shapefile  to  a  grid  and  save  as:  BR_NOGROW. 

3.  Select  limited-access  highways  (Al  class)  from  the  roads  cover  and  make  50m 
buffer  because  developments  right  beside  limited-access  highways  are  regarded 
as  unlikely  convert  to  a  grid:  BR_I_50MBUFF. 

4.  Merge  BRJNTOGROW  and  BR_I_50MBUFF. 

5.  Clip  the  merged  grid  with  the  boundary  grid. 

6.  Reclassify  the  grid  so  that  within  the  area:  1  and  others:  0. 

AML  Script  Opportunity:  Procedures  (2)  thru  (6)  can  be  done  by  running  the  No 
Growth  AML  script:  “Arc:  &r  C:\Fort_Bragg  \  9nogrow.aml” 

-  For  the  input  shape  file  name,  type  “BR_NOGROW.SHP”. 

-  For  the  road  exchange  file  name,  type  “BR_ROAD”. 
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-  For  the  value  item  field,  type  “CLASS”. 

-  For  the  limited-access  highways,  type  “1”. 

-  For  the  input  boundary  cover,  type  “BRJOOUNTY”. 

-  For  the  boundary  grid,  type  “BR  GRID”. 

“Nogrowremap” and  “facility remap”  files  must  be  in  the  workspace 
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8.  Ramp  Attractor  Map.  This  will  describe  taking  a  raw  vector  data  file  of 
roads  in  the  Bragg  area  from  a  data  source,  doing  a  cost  calculation  based  on  the 
travel  time  map  made  previously,  and  manipulating  it  into  a  LEAM  input  map 
file.  The  result  will  be  a  raster  file  showing  travel  time  from  each  grid  cell  to  the 
nearest  highway  ramp  (Figure  22). 


Figure  22.  Travel  time  to  nearest  highway  ramp,  Fort  Bragg  region. 

Source 

Road  network  data  from: 

http://arcdata.esri.com/data/tiger2Q0Q/tiger  download. cfm 
Travel  minute  grid  from  4  (e.g.,  TRAVELMIN ) 

Procedure 

1.  Select  ramps  from  the  roads  cover  and  converted  to  a  grid:  BR_RAMP. 

2.  Run  COSTDISTANCE  command  with  BR_RAMP  as  source  and  TRAVELMIN  as 
cost  and  set  the  values  to  integer. 

AML  Script  Opportunity:  This  entire  procedure  can  be  done  using  the  ramp  at¬ 
tractor  AML  script:1  Arc:  &r  C:\Fort_Bragg\10rampatr.a7nl”. 
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9.  Road  Intersection  Map.  This  will  describe  taking  a  raw  vector  data  file 
of  roads  in  the  Bragg  area  from  a  data  source,  extracting  intersection  nodes  from 
the  file,  and  manipulating  it  into  a  LEAM  input  map  file.  The  result  will  be  a 
raster  file  showing  travel  time  to  major  intersections  (Figure  23). 


Figure  23.  Travel  time  to  major  intersections,  Fort  Bragg  region. 
Source 


Road  network  data  from  #2 

Procedure 

1.  A2  and  A3  classes  are  selected  from  the  road  data  and  saved  as  a  shapefile. 

2.  Generate  a  point  theme  where  the  points  are  intersections  of  the  selected  roads: 
an  avenue  script  was  made  by  Yongwook  Kim. 

3.  Use  the  project  file  intersection_point.apr. 

4.  Add  the  necessary  road  shapefile. 

5.  Press  “C”  icon  to  consolidate  arc  lines  A  a  new  shapefile  is  made. 

6.  Press  “N”  icon  to  generate  a  node  shapefile  from  the  new  shapefile. 

7.  Pay  attention  to  the  pop-up  windows. 

8.  Convert  the  point  theme  to  a  grid:  BR_RO_NODE. 

9.  Run  COSTDISTANCE  command  with  BR_RO_NODE  as  source  and 
TRAVELMIN  as  cost  and  set  the  values  to  integer. 

AML  Script  Opportunity:  Procedures  (3)  thru  (4)  can  be  done  using  the  intersec¬ 
tion  AML  script:  “Arc:  &r  C:  \Fort_Bragg  1 1  lrdintatr.aml” 

Must  have  the  point  shapefile  in  the  workspace 
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10.  Slope  Map.  This  will  describe  taking  a  digital  elevation  model  (DEM) 
of  the  Bragg  area  from  a  data  source,  deriving  the  slope  of  the  land,  and  manipu¬ 
lating  it  into  a  LEAM  input  map  file.  The  result  will  be  a  raster  file  showing 
physical  slope  for  each  grid  cell  in  the  Bragg  area  (Figure  24). 


Figure  24.  Physical  slope,  Fort  Bragg  region. 


Source 

USGS  1:24K  DEM  (the  second  one  is  recommended): 
http://edc.usqs.qov/qeodata/ 

http://seamless.usqs.gov/ 

Procedure 

1.  Go  to  http://seamless.usqs.gov/.  use  download  tool  on  National  Map,  extract  DEM  as  a 
zip  file  ->  Download  file  to  appropriate  folder  (there  may  be  more  than  one  file  to 
download)  If  there  is  more  than  one  DEM  file,  you  will  need  to  merge  them 
into  one  file;  this  can  be  done  in  Arclnfo 

(e.g.,  “Arc:  grid  ->  Grid:  br_demmerge  =  merge(c:\fort_Bragg\XXXXXXXX, 
c:\fort_Bragg\YYYYYYYY)”),  where  XXXXXXXX  and  YYYYYYYY  are  8-digit 
numbers  representing  the  downloaded  DEM  section  files. 

Note:  Due  to  the  download  process,  the  8-digit  DEM  files  may  be  held  in  a  sub¬ 
folder  with  the  same  name  as  the  8-digit  filename.  In  this  case,  the  Arc  com¬ 
mand  would  be  “Arc:  grid  ->  Grid:  br_demmerge  = 
merge(c:\fort_Bragg\XXXXXXXX\XXXXXXXX, 
c:\fort_Bragg\YYYYYYYY\YYYYYYYY)” 

Note:  to  get  out  of  Grid  mode,  enter  “q”  or  “quit” 

2.  Project  merged  grid  to  Albers  projection  (this  will  take  a  long  time  in  Arclnfo). 

3.  Clip  the  grid  with  the  boundary  grid. 
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4.  Derive  slope  using  ArcView  derive  slope  menu. 

AML  Script  Opportunity:  Procedures  (2)  thru  (4)  can  be  done  using  the  slope 
AML  script: 

“Arc:  &r  C:  \Fort_Bragg  \  12slope.aml” 
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11.  State  Road  Attractor  Map.  This  will  describe  taking  a  raw  vector  data 
file  of  roads  in  the  Bragg  area  from  a  data  source,  doing  a  cost  calculation  based 
on  the  travel  time  map  made  previously,  and  manipulating  it  into  a  LEAM  input 
map  file.  The  result  will  be  a  raster  file  showing  travel  time  from  each  grid  cell 
to  the  nearest  state  highway. 


Figure  25.  Travel  time  to  state  highways,  Fort  Bragg  region. 

Source 

Road  network  data  as  an  exchange  file  from  #4 
Travel  minute  grid  from  4  (e.g.,  TRAVELMIN ) 

Procedure 

1.  Select  state  highways  from  the  roads  cover  and  convert  to  a  grid:  BR_STHWYS. 

2.  Run  COSTDISTANCE  command  with  BR_STHWYS  as  source  and 
TRAVELMIN  as  cost  and  set  the  values  to  integer. 

AML  Script  Opportunity:  This  entire  procedure  can  be  done  using  the  state  high¬ 
way  attractor  AML  script:1  Arc:  &r  C:  \Fort_Bragg\  13stateatr.  ami” 

Tree  Attractor  Map.  This  will  describe  taking  information  from  the  Bragg  area 
NLCD  map,  calculating  distance  to  the  nearest  forest,  and  manipulating  it  into  a 
LEAM  input  map  file.  The  result  will  be  a  raster  file  showing  Euclidean  dis¬ 
tance  from  each  grid  cell  to  the  nearest  forest  grid  cell  (Figure  26). 
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Figure  26.  Distance  to  nearest  forest,  Fort  Bragg  region. 

Source 


USGS  NLCD 

Procedure 

1.  Select  forest  cover  cells  (41,  42  and  43)  from  the  NLCD  grid:  BR_FOREST. 

2.  Use  FIND  DISTANCE  menu  to  find  distance  to  forest  cells  and  set  the  values  to 
integer. 

AML  Script  Opportunity:  This  entire  procedure  can  he  done  using  the  tree  attrac¬ 
tor  AML  script: 

“Arc:  &r  C:  \Fort_Bragg  \14treeatr.aml”. 
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12.  Water  Attractor  Map.  This  will  describe  taking  information  from 
the  Bragg  area  NLCD  map,  calculating  distance  to  the  nearest  water  body,  and 
manipulating  it  into  a  LEAM  input  map  file.  The  result  will  be  a  raster  file 
showing  Euclidean  distance  from  each  grid  cell  to  the  nearest  water  or  wetland 
grid  cell  (Figure  27). 


Figure  27.  Distance  to  nearest  water,  Fort  Bragg  region. 

Source 


USGS  NLCD 

Procedure 

1.  Select  open  water  and  wetland  cells  (11,  91,  and  92)  from  the  NLCD  grid: 
BRWATER. 

2.  Use  FIND  DISTANCE  menu  to  find  distance  to  water  cells  and  set  the  values  to 
integer. 

AML  Script  Opportunity:  This  entire  procedure  can  he  done  using  the  water  at¬ 
tractor  AML  script:  “Arc:  &r  C:\Fort_Bragg  1 15wateratr .ami” 

•k’k'k'k'k'k'k’k'k'k'k’k'k'k'k'k’k'k'k'k'k’k'k'k'k'k'k'k'k'k'k’k'k’k'k'k'k'k'k'k’k'k'k'k'k’k'k’k'k'k'k'k'k'k'k'k'k'k'k’k'k'k'k'k’k 

All  maps  created  in  this  tutorial  are  mLEAM  input  maps  and  are  now 
ready  for  use  with  the  LEAM  software. 

•k’k’k’k’k'k'kic'k’k'kirk’k’k'k’k’k’k’kic'k’k’k'k'k’k’k’k’k’k'k’k’k'k’k’k’kic’k’k’k’kic'kic’k’k’k'kic'k'k'k'kic’k’k’k’kic'kic'k'k 
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ESRI  ArcGIS  with  GRASS  Version 

Initial  Setup 

1.  On  Windows  machine,  create  directory  on  hard  drive  (e.g.,  C:\Fort_Bragg)  to 
hold  data.  There  should  not  be  any  spaces  in  the  path  names. 

2.  On  Windows  machine,  install  the  program  WinSCP — this  will  allow  the  user  to 
move  files  back  and  forth  between  the  Windows  environment  and  the  UNIX  en¬ 
vironment.  WinSCP  is  available  for  download  on  the  internet. 

3.  On  Windows  machine,  ensure  that  ESRI  ArcGIS  is  installed. 

4.  On  UNIX  computer  (or  an  X-Windows  interface  to  a  remote  UNIX  computer) , 
ensure  that  GRASS  is  installed. 

5.  Before  the  manipulated  data  is  brought  into  GRASS,  it  is  necessary  to  create  two 
different  LOCATIONS  in  GRASS.  One  LOCATION  should  have  a  Latitude- 
Longitude  Geographic  coordinate  system  and  should  be  named  as  such.  The 
other  LOCATION  should  have  an  Albers  Equal  Area  projection  system  and 
should  be  named  as  such.  These  two  LOCATIONS  will  make  it  possible  to  im¬ 
port  most  of  the  manipulated  data  files,  as  these  files  typically  are  projected  in  ei¬ 
ther  Lat/Long  or  Albers.  The  r.proj  command  can  take  a  file  from  one 
LOCATION  and  reproject  it  into  the  current  LOCATION’S  projection. 

Creating  Base  Maps 

Mapl:  County  boundary  map  (BR  COUNT).  The  GRASS  version  is  a  raster  file 

showing  the  counties  of  interest. 

Creating  Map  1  from  scratch  using  ArcToolbox  and  ArcMap  8.3: 

Source:  County  boundary  data  from  URL: 

http://www.census.qov/qeo/www/cob/co2000.html 

1.  Select  North  Carolina,  download  and  unzip  the  .eOO  file  to  the  directory  of  choice. 

2.  Import  .eOO  file  using  ArcToolbox:  Conversion  Tools -^Import  to  cover¬ 
age  -^Import  from  interchange  file. 

Input  file=<co37_dOO.eOO> 

Output  dataset=<br_county>  (Ensure  the  output  file  path  is  the  same  as 
the  input  file  path.) 

3.  Project  file  to  Albers  using  Projections-^  Project  wizard. 

4.  Convert  to  shapefile  using  ArcToolbox:  Conversion  Tools-> Export  from  Cover¬ 
age-^  Coverage  to  Shapefile. 

5.  Select  necessary  counties  (Cumberland,  Harnett,  Hoke,  Lee,  Moore,  Richmond, 
Scotland)  using  ArcMap,  then  save  the  selection  as  a  new  file. 
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6.  Copy  this  file  to  UNIX  environment  using  WinSCP. 

7.  The  resulting  map  should  be  a  shapefile  of  the  counties  of  interest,  and  this  can 
be  imported  into  GRASS  using  r.in.shape.  The  GRASS  location  should  be  the 
one  with  the  Albers  Equal  Area  projection  defined. 

(Example:  r.in.shape  in=counties.shp  out=counties) 

Map  2:  NLCD  land  use  data  (BR  NLCD).  Creating  Map  2  from  scratch: 

Source:  USGS  National  Land  Cover  Data  (NLCD)  from  URL: 
http://seamless.usqs.qov 

1.  Download  a  NLCD  TIFF  file  of  North  Carolina  and  unzip. 

2.  Save  as  grid  (nc_nlcd.tif)  and  copy  to  UNIX  environment  using  WinSCP. 

3.  This  file  can  be  imported  into  GRASS  using  r.in.gdal. 

(Example:  r.in.gdal  — o  input=nc_nlcd.tif  output=nc_nlcd) 

Map  3:  Roads  (BR_ROADS).  GRASS  version  is  raster  file  showing  roads  on 
counties  of  interest.  The  roads  are  classified  according  to  speed  limits  (Al 
through  A7). 

Creating  Map  3  from  scratch: 

Source:  Road  network  data  from  Census: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 

1.  Download  road  shapefiles  of  counties  of  interest  (Cumberland,  Harnett,  Hoke, 
Lee,  Moore,  Richmond,  Scotland)  and  merge  into  one  file  in  ArcMap  8.3  (e.g., 
ROAD.SHP). 

2.  Open  ArcView  3.2  and  make  a  new  field  for  the  roads: 

Open  ROADS  theme  and  open  the  attribute  table. 

Create  two  new  fields,  one  titled  “Class”  and  one  titled  “New.”  “Class” 
should  be  a  numeric  field,  and  “New”  should  be  a  text  field. 

Start  editing. 

Select  all  entries  in  the  attribute  table. 

Calculate  values  for  “New”  with  the  expression  =[Cfcc2].Right(l). 

The  switch  for  String  will  have  to  be  checked. 

Calculate  values  for  “Class”  with  the  expression  =  [Ne wj . As N u  m b e r . 

The  result  should  be  that  “Class”  shows  the  numerical  values  for  “Cfcc2,” 
without  the  “A”  character  in  front  of  the  number. 

Save  edits. 

3.  Copy  this  file  to  UNIX  environment  using  WinSCP. 
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4.  This  file  can  be  imported  into  GRASS  using  r.in. shape.  It  is  important  to  include 
the  additional  parameter  “cat=class”  so  that  the  roads  are  classified  according  to 
their  speed  limit  classes.  (Example:  r.in.shape  in=road.shp  out=road  cat=class). 


Map  4:  Municipal  boundaries  map  (BR  MUNBND).  GRASS  version  is  raster  file 
showing  municipal  boundaries  of  all  towns  within  the  counties  of  interest. 

Creating  Map  4  from  scratch: 

Source: 

Municipal  boundary  data  from  Census 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 

Census  population  data  of  incorporated  places:  Summary  File  l->Geographic 
Comparison  Tables-^ State->Place  (GCT-PH1) 

1.  Download  municipal  boundary  shapefiles  of  counties  of  interest  (Cumberland, 
Harnett,  Hoke,  Lee,  Moore,  Richmond,  Scotland)  and  merge  into  one  file  in 
ArcMap  8.3  (e.g.,  MUNBND.SHP). 

2.  Open  the  table  of  the  shapefile,  add  a  field  “population”  with  long  integer,  start 
editing,  and  input  the  values  from  Census.  Save  edits. 

3.  Copy  this  file  to  UNIX  environment  using  WinSCP. 

4.  This  file  can  be  imported  into  GRASS  using  r.in.shape.  To  create  maps  5,  6,  7,  it 
is  also  important  to  include  the  additional  parameter  “cat=Population.”  This  will 
probably  need  to  be  imported  into  the  Albers  LOCATION,  since  the  merged  mu¬ 
nicipal  boundary  layers  usually  are  projected  with  Albers  by  default. 

(Example:  r.in.shape  in=munbnd.shp  out=munbnd  cat=  Population) 

Maps  5,  6,  7.  (Small  city  boundaries  map  (BR  SMCITY);  Medium  city  boundaries 
map  (BR  MEDCITY);  Large  city  boundaries  map  (BR_LGCITY).  GRASS  version  is 
three  distinct  raster  files  showing  municipal  boundaries  of  small,  medium,  and 
large  cities  within  the  counties  of  interest.  These  files  were  created  from  the 
BR_MUNBND  map  (map  4)  above.  The  cities  were  arbitrarily  classified  into 
three  groups  with  population  data  from  the  U.S.  Census. 

Creating  Maps  5,  6,  7  from  scratch: 

If  the  BR_MUNBND  map  (Map  4)  was  imported  with  the  category  “Population,” 
maps  5,  6,  and  7  can  be  created  in  GRASS. 

1.  In  GRASS,  use  r.reclass  to  categorize  the  city  polygons  into  three  classes  based 
on  their  population.  The  result  will  be  a  new  GRASS  map. 

(Example:  r.reclass  munbnd  out=class_munbnd 
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1  thru  1000  =  1  small  city 
1001  thru  24000  =  2  med  city 
24001  thru  150000  =  3  large  city) 

2.  In  GRASS,  use  r.reclass  again  (this  time,  on  the  new  map  created  in  step  2)  to 
make  a  map  of  small  city  polygons  only. 

(Example:  r.reclass  class_munbnd  out=small_city 

1  =  1  small  city 

2  thru  3  =  NULL) 

3.  Repeat  step  3  to  make  maps  for  medium  city  polygons  and  large  city  polygons. 

4.  The  result  will  be  three  different  maps,  each  showing  polygons  of  small,  medium, 
or  large  cities  (maps  5,  6,  7). 

Map  8.  No  growth  map  (BR  NOGROWTH).  GRASS  version  is  raster  file  show¬ 
ing  areas  where  growth  will  not  occur.  These  areas  include  U.S.  military  instal¬ 
lations,  Government  lands,  landmarks,  water  bodies,  Indian  lands,  and  the  like. 

Creating  Map  8  from  scratch: 

Source  (use  some  or  all  of  these  sources) 

National  Atlas  for  Federal  and  Indian  land: 
http://www.nationalatlas.gov 

Other  local  data  about  ownership 

Landmark  polygon  data  from  URL: 

http://arcdata.esri.com/data/tiger200Q/tiger  download. cfm 

ROADS  map  (Map  #3) 

Boundary  cover  shapefile 

1.  In  ArcMap  8.3,  gather  the  appropriate  maps  for  the  counties  of  interest  from 
sources  i  through  iii  above  and  merge  into  one  “no  growth”  shapefile.  You  will 
need  to  clip  the  Federal/Indian  map  with  the  county  boundary  shapefile  before 
merging  with  the  other  maps. 

2.  Copy  this  file  to  UNIX  environment  using  WinSCP. 

3.  This  file  can  be  imported  into  GRASS  using  r.in.shape. 

4.  (Example:  r.in.shape  in=no_grow  out=no_grow). 

5.  In  GRASS,  use  r.reclass  to  make  a  new  road  map  with  just  class  A1  roads. 

(Example:  r.reclass  roads  out=interstate 

1  =  1  interstate 

2  thru  7  =  NULL) 
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6.  Use  r. buffer  to  create  a  50m  buffer  around  the  Al  roads,  because  growth  closer 
than  50m  to  an  interstate  is  highly  unlikely. 

(Example:  r. buffer  interstate  out=interstate_buffer  dist=50  unit=meters) 

7.  In  GRASS,  merge  the  no  growth  shapefile  and  the  interstate  buffer  file  using 
r.mapcalc. 

(Example:  r.mapcalc 

‘nogrowth=if(isnull(interstate_buffer),if(no_grow),if(interstate_buffer))’ 

8.  This  will  make  a  new  map  showing  areas  of  EITHER  50m  interstate  buffer  OR 
Federal/Indian/water  lands. 

Map  9.  Digital  elevation  model  (BR  DEM)  GRASS  version  is  raster  file  showing 
elevation  of  the  area  of  interest.  The  map  will  contain  the  counties  of  interest 
but  will  also  have  extraneous  land  because  of  the  way  the  DEM  data  was  ex¬ 
tracted  (the  file  will  not  have  the  “cookie  cutter”  shape  of  the  county  boundary 
map,  it  will  be  a  rectangular  grid  containing  the  counties  of  interest). 

Creating  Map  9  from  scratch: 

Source:  USGS  1:24K  DEM.  The  second  URL  (the  National  Map)  is  recom¬ 
mended: 

http://edc.usqs.qov/qeodata/ 

http://seamless.usqs.qov/ 

1.  To  obtain  the  DEM  from  The  National  Map,. 

Go  to  http://seamless.usqs.qov/ ,  zoom  in  to  the  area  of  interest 

Click  on  “NED”  checkbox  in  upper  right  hand  menu  under  Download 
Layers.  Uncheck  all  other  boxes. 

Use  rectangular  download  tool  (the  icon  is  on  the  middle  left  hand  side) 
and  draw  around  the  area  of  interest. 

The  website  should  automatically  begin  downloading  and  might  take 
awhile  to  finish. 

2.  Unzip  the  downloaded  file  and  save  to  the  appropriate  directory. 

3.  Copy  the  finished  grid  file  to  UNIX  environment  using  WinSCP. 

4.  This  grid  can  be  imported  into  GRASS  using  r.in.gdal  or  r.in.arc. 

(Example:  r.in.gdal  in=dem  out=dem) 

Warning:  the  grid  file  may  have  an  eight  digit  name  assigned  to  it  due  to 
the  National  Map  conventions. 
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Creating  mLEAM  Input  Maps  via  GRASS  Makefile: 

From  the  nine  base  maps  created  above,  the  GRASS  makefile  can  automatically 
produce  the  mLEAM  input  maps.  The  code  for  the  GRASS  makefile  is  shown  in 
Appendix  C. 


Table  2.  mLEAM  input  files  made  by  GRASS  makefile. 


Map 

Description 

Boundary 

Boundary  of  study  region  (usually  a  group  of  counties) 

Land  Cover 

Land  cover  of  study  region  (NLCD  data) 

Metrobuffer 

Buffer  around  municipal  boundaries  in  study  region  (typ  1 .5  miles) 

No  Growth 

Areas  in  study  region  in  which  development  is  not  possible  or  highly 
unlikely 

Slope 

Slope  of  terrain  in  study  region 

Forest  Attractor 

Distance  to  nearest  forest 

Water  Attractor 

Distance  to  nearest  water  body 

Highway  Attractor 

Travel  time  to  nearest  highway 

Intersection  Attractor 

Travel  time  to  nearest  major  road  intersection 

Ramp  Attractor 

Travel  time  to  nearest  highway  ramp 

Road  Attractor 

Travel  time  to  nearest  road 

Cities  Attractor  (large, 
medium,  small) 

Travel  time  to  nearest  city  (divided  into  three  groups — small,  medium, 
and  large) 
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Create  Alternative  Scenarios 

A  new  loop  highway,  running  from  NW  Fayetteville  along  the  edge  of  Fort  Bragg 
and  connecting  with  1-95  again  in  SW  Fayetteville,  was  scheduled  to  be  built  in 
the  near  future.  The  CAD  file  for  this  highway  was  provided  by  Fort  Bragg,  and 
then  digitized  into  GIS  using  ArcMap  8.3.  This  road  was  then  added  to  the 
original  road  network  map,  which  resulted  in  a  secondary  road  network  map; 
this  acted  as  a  base  map  for  a  different  mLEAM  growth  scenario.  Two  different 
possibilities  for  the  area  emerged:  (1)  without  the  new  highway,  or  (2)  with  the 
new  highway. 


Table  3.  Additional  mLEAM  input  files  made  by  GRASS  makefile  due  to  new  highway. 


Map 

Description 

Cities  Attractor  (large, 
medium,  small) 

Travel  time  to  nearest  city  (divided  into  three  groups — small,  medium, 
and  large) 

Highway  Attractor  2 

Travel  time  to  nearest  highway 

Intersection  Attractor  2 

Travel  time  to  nearest  major  road  intersection 

Ramp  Attractor  2 

Travel  time  to  nearest  highway  ramp 

Road  Attractor  2 

Travel  time  to  nearest  road 
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Run  mLEAM 

The  resulting  mLEAM  input  maps  can  be  put  into  the  LEAM  software  to  pro¬ 
duce  the  final  mLEAM  output  models.  The  results  of  running  this  model  are 
shown  in  Figures  28  to  31.  The  two  resultant  maps  showed  growth  for  the  two 
different  road  scenarios:  (1)  without  the  new  highway,  and  (2)  with  the  new 
highway. 


Figure  28.  futureLanduse@scenarioBase. 
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Figure  29.  futureLanduse@scenarioRoad. 


Figure  30.  ATTRACTOR_RES@scenarioBase 
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Figure  31.  ATTRACTOR_RES@scenarioRoad. 


Review  and  Compare  Results 

Each  of  the  mLEAM  outputs  showed  considerable  growth  of  Fayetteville  toward 
Fort  Bragg,  with  development  touching  the  border  of  the  installation  in  many 
cases.  Large  growth  was  also  predicted  on  the  western  edge  of  Fort  Bragg,  cen¬ 
tered  on  the  cities  of  Pinehurst  and  Southern  Pines. 

The  growth  patterns  around  the  Fayetteville  region,  however,  differed  between 
the  two  scenarios.  The  scenario  without  the  new  highway  showed  a  more  cen¬ 
tralized  growth  pattern  situated  primarily  along  the  main  arteries  going  in  and 
out  of  Fayetteville.  The  scenario  with  the  new  highway  predicted  more  growth 
along  this  highway,  resulting  in  much  more  growth  closer  to  Fort  Bragg  and  to 
the  north  of  Fayetteville.  This  is  expected,  since  residents  usually  desire  to  live 
close  to  roads  that  have  quick  access  to  and  from  their  places  of  work.  A  map 
was  made  (Figure  32)  highlighting  the  difference  in  the  growth  patterns  between 
the  two  scenarios  (this  map  was  not  made  in  the  mLEAM  program — it  was  done 
manually  in  a  GIS).  Figures  33  to  36  show  the  schematics  for  converting  the 
map  data. 
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Figure  32.  Difference  in  the  growth  patterns  between  two  scenarios. 


Summary  of  Data  Sources 

NC  County  boundary  data: 

http://www.census.gov/qeo/www/cob/co2000.html 

USGS  NLCD  data  for  NC: 

http://edcsqs9.cr.usqs.gov/pub/data/landcover/states 

NC  roads  by  county: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 
NC  municipal  boundaries: 

http://arcdata.esri.com/data/tiqer2000/tiqer  download. cfm 

Census  population  NC  data: 

http://www.census.gov 

NC  Federal  &  Indian  lands: 

http://www.nationalatlas.gov 

NC  landmark  polygons: 

http://arcdata.esri.com/data/tiger2000/tiqer  download. cfm 

NC  environmental  lands:  JLUS  CD  (provided  by  Fort  Bragg  area  commission) 
USGS  1:24K  DEM: 

http://seamless.usqs.qov 

New  highway:  provided  by  Fort  Bragg,  digitized  by  Kyle  Brock 
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N.  C.  county  polygon 
(from  Census) 

co37_d00.e00 

> 


> 


Fort  Bragg  area  counties 
(coverage) 
BFLCOUNTY 


Fort  Bragg  area  counties 

(grid) 

BFI.GRID 


These  2  files  are  also 
used  as  intermediates  for 
making  other  LEAM  files. 


North  Carolina  NLCD 

Fort  Bragg  NLCD  grid 

(from  USGS) 

BR_NLCD 

North  Carolina  roads  by 
county 
(from  ESRI) 


Fort  Bragg  area  roads 
BR_ROAD 


Ft  Bragg  area  tree 
attractor 
BR_TREEATR 


Ft  Bragg  area  water 
attractor 

BR  WATERATR 


Fort  Bragg  land  use 
BR_NLCD_R 


Travel  time  between 
cells,  Ft  Bragg  area 
BR  TRAV 


'Ft  Bragg  area  major 
road  intersections 
BR  INT 


Ft  Bragg  area  road 
attractor 
BR_ROADATR 


Ft  Bragg  area  small  city 
attractor 
BR_RAMPATR 


North  Carolina  population 
data 

(from  Census) 

N.  C.  municipal  bounds 
by  county 
(from  ESRI) 

Ft  Bragg  area  muni  bnds 
of  small  cities  BR_SCITY 


Ft  Bragg  area  state 
highway  attractor 
BR_STRDATR 


Ft  Bragg  area  intersect, 
attractor 
BRJNTATR 


Ft  Bragg  area  muni  bnds 
of  med  cities 
BR_MCITY 


Ft  Bragg  area  small  city 
attractor 
BR_SMCIT 


Ft  Bragg  area  muni  bnds 
of  large  cities 
BR_LCITY 


Ft  Bragg  area  medium 
city  attractor 
BR.MDCIT 


Other  local  data 
(any  other  sources) 

Federal/Indian  lands 
polygons 
(from  Nat’l  Atlas) 

"i 

Landmark  polygons 
(from  ESRI) 

Ft  Bragg  area  municipal 
boundaries 
BR_MUNBND 


Ft  Bragg  no  growth  areas 
BR_NOGROW 


Ft  Bragg  area  large  city 
attractor 
BR_LGCIT 


Ft  Bragg  municipal 
boundary  buffer 
BR_MUNBUFF 


Ft  Bragg  no  growth  areas 
BR_NOGROWTH 


Dotted  lines  denote  processes 
done  by  hand. 

Solid  lines  denote  processes 
done  through  automatic 
processing  (AML  scripts). 


Ft  Bragg  area  DEM 
(from  USGS) 
BR_DEM 


Ft  Bragg  area  terrain 
slope 

BR_SLOPE 


:igure  33.  Schematic  of  data  conversion  using  ArcGIS. 
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Figure  34.  Schematic  of  data  conversion  using  Grass. 
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Figure  35.  Schematic  of  data  conversion  using  Arc/Info  only  (new  highway  scenario). 
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Figure  36.  Schematic  of  data  conversion  using  ArcGIS  and  Grass  (new  highway  scenario). 
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5  Conclusion 


Prediction  of  urban  patterns  is  accomplished  using  the  Land  Use  Change  (LUC) 
model,  which  is  part  of  a  suite  of  tools  called  the  military  Land  use  Evolution 
and  Assessment  Model  (mLEAM).  This  work  has  provided  information  to  assist 
GIS  technicians  in  preparing  input  maps  for  the  mLEAM  LUC  model,  and  also 
descriptions  of  the  specific  steps  used  to  create  maps  for  the  multi-county  area 
surrounding  Fort  Bragg,  NC.  This  provides  installation  planners  with  better 
tools  to: 

1.  Predict  urban  growth  patterns  that  might  occur  in  response  to  proposed  regional 
plans  involving  certain  investments  and  policies,  and 

2.  Predict  the  change  in  training  and  testing  opportunity  that  might  occur  as  a  re¬ 
sult  of  those  patterns. 
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Appendix  A:  CD  Contents  and  Directory 

Listing 

CD — Fort  Bragg 

The  CD  for  Fort  Bragg  mLEAM  includes  all  files  associated  with  the  production 
of  the  mLEAM  scenario.  Table  Al  lists  the  files.  (You  will  need  to  unzip  (un¬ 
compress)  the  file  on  the  CD  to  read  the  files. 

File  Folders  (The  files  in  these  folders  are  largely  support  or  starter  files  for  the 
maps  shown  in  the  table  below.  There  is  little  reason  to  look  within  these  fold¬ 
ers): 

•  calc 

•  Clip 

•  Tiger_files 

•  Info. 


Table  Al.  Files  associated  with  the  production  of  the  mLEAM  scenario. 


File  name 

Description 

How  obtained 

br_county 

Polygon  coverage  of  7-county  area  around  Fort  Bragg. 
Includes  Cumberland,  Harnett,  Hoke,  Lee,  Moore, 
Richmond,  and  Scotland  Counties. 

AML  script 
“1bnd_county.aml” 

co37_d00 

North  Carolina  county  boundaries,  projected  Albers 

internet 

55687211 

Bragg  area  DEM — Western  half 

internet 

87615094 

Bragg  area  DEM — Eastern  half 

internet 

92088612 

Bragg  area  DEM — entire 

internet 

nc_nlcd 

NLCD  for  North  Carolina 

manipulated  from  Tiff 

north_carolina_ 

NLCD_erd04220.tif 

Tiff  image  of  NLCD  of  North  Carolina 

internet 

br_dem 

Bragg  area  DEM 

manipulated  from  8- 
digit  DEM 

br_demjunk2 

Bragg  area  DEM,  clipped  with  counties 

manipulated  from 
br_dem  and  br_county 

br_demmerge 

Bragg  area  DEM,  entire 

combination  of  W  and  E 

halves 
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File  name 

Description 

How  obtained 

br_nlcd 

NLCD  data  for  Fort  Bragg  area 

AML  script 
“Ibnd  grid. ami” 

brgrid 

Raster  grid  of  7-county  area 

“Ibnd  grid. ami” 

brju 

NLCD  data  for  Fort  Bragg  area,  with  major  roads 
merged 

“2landuse_new.aml” 

br_buf 

7-county  area — 1 .5  mi  buffer  around  municipal  bounda¬ 
ries 

“3munbuf.aml” 

br_trav 

Travel  time  across  grid  cell;  “travel  time”  is  based  on  the 

“4  travel,  a  ml” 

br_roadatr 

brjgcit 

br_mdit 

br_smcit 

br_nogrowth 

br_rampatr 

brjntatr 

br_slope 

br_stateatr 


br_treeatr 
br  wateratr 


number  of  minutes  it  takes  to  travel  across  one  grid  cell, 
and  is  based  solely  on  that  cell’s  NLCD  land  cover  cate¬ 
gory  and/or  existence  of  roads 

Travel  time  to  nearest  road _ “6roadatr.aml” 

Travel  time  to  nearest  large  city  “8cityatr.aml” 

Travel  time  to  nearest  medium  sized  city  “8cityatr.aml” 

Travel  time  to  nearest  small  city _ “8cityatr.aml” 

No  growth  areas  in  Fort  Bragg  area  “9nogrow.aml” 

Travel  time  to  nearest  highway  ramp _ “IQrampatr.aml” 

Travel  time  to  nearest  major  intersection  “1  Irdintatr.aml” 

Physical  slope  of  each  grid  cell  “12slope.aml” 

Travel  time  to  nearest  state  highway  “13stateatr.aml” 


Distance  to  nearest  forest  “14treeatr.aml” 

Distance  to  nearest  water/wetland  “15wateratr.aml” 


a2a3road.shp 

brjcity.shp 


br_mcity.shp 

br_scity.shp 

br_munbnd.shp 

nogrowth.shp 

br_road.shp 

Bragg_bnd.shp 

Bragg_county.shp 

brng.shp 

cn_a2a1  .shp 

nd_a2a1.shp 


Major  roads,  used  in  determining  intersections 

Boundaries  of  large  cities  within  Fort  Bragg  area — 
shapefile 


Boundaries  of  medium  sized  cities  within  Fort  Bragg 
area — shapefile 

Boundaries  of  small  cities  within  Fort  Bragg  area — 

shapefile _ 

Boundaries  of  cities  within  Fort  Bragg  area 

No  growth  areas  in  Fort  Bragg  area 

All  roads  in  Fort  Bragg  area 

Boundary  of  Fort  Bragg 

7-county  area,  projected  Albers 

No  growth  areas  in  Fort  Bragg  area 

Major  roads,  used  in  determining  intersections 

Major  road  intersections 


manipulated  from 
br_road.shp  using  Arc- 
View  3.2 

manipulated  from 
br_munbnd.shp 


manipulated  from 
br_munbnd.shp 
manipulated  from 

br_munbnd.shp _ 

combined  from  Census 
Tiger  files 

combined  from  various 
sources 

combined  from  Census 
Tiger  files 

selected  from  U.S.  mili¬ 
tary  base  shapefile 
manipulated  from 
co37_d00 

combined  from  various 
sources 

calculated  inArcView 

3.2 _ 

calculated  in  ArcView 
3.2 
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File  name 


nd_cn_2.shp 


Description 


Major  road  nodes  (intersections,  bends) 


How  obtained 


calculated  in  ArcView 
3.2 


fedlanp020.shp 

nc_mil.shp 

nc_state.shp 

Igcityl.shp 

managed. shp 


Federal  &  Indian  lands  internet 

Fort  Bragg  and  Camp  Mackall  boundaries  manipulated  from  man- 

_ aged. shp _ 

N.  Carolina  state  outline  polygon  internet 

Boundaries  of  major  cities  in  Fort  Bragg  area  (similar  but  manipulated  from 
not  the  same  as  brjcity.shp)  br_munbnd.shp 

Managed  environmental  lands  within  7-county  area  combined  from  files  on 

a  CD  donated  by  Fort 
Bragg  planning  com- 


br_count.e00 
br_road1  .eOO 
br_roadnew.eOO 
co37  dOO.eOO 


7-county  area — interchange  file _ 

Bragg  area  roads,  not  projected — interchange  file 
Bragg  area  roads,  projected  Albers — interchange  file 
NC  counties — interchange  file 


“Ibnd  cover.aml” 


internet 


AML  scripts 


1bnd_cover.aml 

1bnd_grid.aml 

2landuse_new.aml 

3munbuf.aml _ 

4travel.aml 


6roadatr.aml 

8cityatr.aml _ 

9nogrow.aml 

lOrampatr.aml 

llrdintatr.aml 


12slope.aml 

13stateatr.aml 


14treeatr.aml 

15wateratr.aml 


Support  files  for  AML 
scripts 
boundremap 
boundremap2 

cityremap _ 

facilityremap 

fid  remap _ 

growremap 


munbndremap 

nogrowremap 


road  remap 
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Appendix  B:  AML  Scripts 

lbnd_cover . ami 

J  ' ************************************************************* 

/*  lbnd_cover . ami  * 

/  -k  -k 

/*  This  AML  will  import  the  *.e00  file  as  cover  and  * 

/*  select  counties  of  interest  with  appropriate  projection  * 

/  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Dec  16,  2002,  * 

/*  modified  by  Woonsup  Choi  on  June  24,  2003  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/***************************** 

/*  setup  the  input  parameter  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

S setvar  input_exchangein  =  [response  'Enter  a  exchange  file  name  without 
extension ' ] 

Ssetvar  cty_bnd  =  [response  'Enter  a  name  of  the  resulting  cover'] 

Ssetvar  input_exchange  =  [locase  %input_exchangein% ] 

&TYPE 

&Type  All  THE  EXISTING  COVERAGE  OR  GRID  NAMED  " %input_exchange% , " 
"%input_exchange%l , "  "%input_exchange%2 , "  "%cty_bnd%"  WILL  BE  DELETED  !!! 
&TYPE 

Ssetvar  delete_yesnoin  =  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

Sif  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  sthen 
Sreturn  ERROR!  You  should  type  either  'y'  or  'n' 

Sif  %delete_yesno%  =  n  sthen 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

/*  Sif  A  [exists  %tifname%  -file]  Sthen 

/*  Sif  [exists  [value  temppr j %pr j % ]  -file]  Sthen 

/‘Ssetvar  input_unitin :=[ response  'Is  the  unit  of  input  file  feet  or  meter? 
(f/m) ' ] 

/‘Ssetvar  input  unit  =  [locase  %input_unitin% ] 

/*Sif  %input_unit%  ne:=f  and  %input_unit%  ne:=m  Sthen 

/*  Sreturn  ERROR!  You  have  to  type  either  'f'  of  'm'  for  the  unit. 

Sif  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

kill  %input_exchange%  all 

Sif  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

Sif  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 

Sif  [exists  %cty_bnd%  -COVER]  or  [exists  %cty_bnd%  -GRID]  Sthen 
kill  %cty_bnd%  all 

Sif  [exists  %cty_bnd%l  -COVER]  or  [exists  %cty_bnd%l  -GRID]  Sthen 
kill  %cty_bnd%l  all 
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/******************************************** 

/*  import  Arclnfo  exchange  file  to  Coverage  * 
/******************************************** 
import  cover  %input_exchange% . eOO  %input_exchange% 

/**************************************************** 

/*  projecting  the  Coverage  to  the  target  projection  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

project  cover  %input  exchange%  %input  exchange%l 
output 

projection  albers 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 
29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 
0 

end 

clean  %input_exchange%l 

j  ' k-k-k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k'k'k-k-k-k-k-k-k'k'k-k-k'k'k-k-k-k'k-k'k-k-k-k-k-k'k-k-k-k 

/*  renaming  the  cover  * 

/*********************************************** 
kill  %input_exchange% 

rename  %input  exchange%l  %input  exchange% 

/********************************************* 

/*  selecting  counties  of  interest  * 

/  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

reselect  %input_exchange%  %cty_bnd%  poly 
lbnd_gr id . ami 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  lbnd_grid. ami  * 

/  -k  -k 

/*  This  AML  will  change  the  coverage  to  GRID  * 

/*  * 

/*  Made  by  Yong  Wook  Kim  on  Dec  16,  2002,  * 

/*  modified  by  Woonsup  Choi  on  July  03,  2003  * 

j  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 
j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  setup  the  input  parameter  * 

/***************************** 

Ssetvar  input_exchangein  =  [response  'Enter  a  county  boundary  cover  name'] 
Ssetvar  input_snlcd  =  [response  'Enter  a  state  NLCD  grid  name  you  have'] 

/*  You  must  have  a  state  NLCD  grid  converted  from  TIFF  file 

Ssetvar  input_bnlcd  =  [response  'Enter  a  boundary  NLCD  grid  name  you  want'] 
/*  This  is  NLCD  of  the  study  area.  This  will  be  used  to  define  boundary  grid 
cell  alignment. 

Ssetvar  cty_bnd  =  [response  'Enter  a  name  of  the  resulting  grid'] 

Ssetvar  input_exchange  =  [locase  %input_exchangein% ] 

STYPE 

SType  All  THE  EXISTING  COVERAGE  OR  GRID  NAMED  " %input_exchange% 1 , " 
"%input_exchange%2 , "  "%input_bnlcd% , "  "%cty_bnd%, "  "%cty_bnd% "_1  WILL  BE 
DELETED  ! ! ! 

STYPE 
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Ssetvar  delete_yesnoin  =  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 
Sreturn  ERROR!  You  should  type  either  ' y'  or  ' n' 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

&if  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 

&if  [exists  %input_bnlcd%  -COVER]  or  [exists  %input_bnlcd%  -GRID]  Sthen 
kill  %input_bnlcd%  all 

&if  [exists  %cty_bnd%  -COVER]  or  [exists  %cty_bnd%  -GRID]  Sthen 
kill  %cty_bnd%  all 

&if  [exists  %cty_bnd%l  -COVER]  or  [exists  %cty_bnd%l  -GRID]  Sthen 
kill  %cty_bnd%l  all 

&if  [exists  %cty_bnd%_l  -COVER]  or  [exists  %cty_bnd%_l  -GRID]  &then 
kill  %cty_bnd%_l  all 

/*polygrid  %input_exchange%  %cty_bnd%l 
/*30 

/*y 

grid 

gridclip  %input_snlcd%  %input_bnlcd%  cover  %input  exchange% 

%cty_bnd%  =  %input_bnlcd%  *0+1 

q 

gridascii  %cty_bnd%  %cty_bnd% . asc 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  exporting  the  cover  as  eOO  * 

/*********************************************** 
export  cover  %input  exchange%  %input  exchange% . eOO 

21anduse_new . ami 

/* - 


/*  YWKIM's  AML  for  Fort  Bragg 
/* - 


/*  Program:  21anduse_new . ami 

/*  Purpose:  Convert  the  ESRI  street  cover  to  GRID  including  A1-A3,  and  A6. 
/* - 


/*  Usage:  21anduse_new  <INSTR_SHP>  <INLU_GRID>  <INMUN_SHP>  <OUT_GRID> 

/* 

/*  Arguments:  <INSTR_SHP>  -  INPUT  SHAPEFILE  THAT  CLIPPPED  FROM  ESRI  STREET 
MAP 

/*  <INLU_GRID>  -  INPUT  NLCD  LANDUSE  GRID  of  the  study  area  —  should 

be  clipped  from  the  state  data 
/*  <OUT_GRID>  -  OUTPUT  GRID 

/*  <INMUN_SHP>  -  INPUT  MUNICIPAL  BOUNDARY  FROM  TIGER 

/* - 


/* 

/* 

/* 


Prerequisite:  Street  Coverage  converted  from  the  ESRI  street  map 
Process:  Clip  the  shapefile,  shapearc,  then  routearc 


/************************* 
/*  basic  parameter  setup  * 


ERDC/CERL  TR-04-24 


65 


/************************* 

&args  incover  lugrid  inmun  outgrid 

&if  [NULL  %incover%]  =  .TRUE.  &then 

Sreturn  Sinform  Usage:  ASTRGRD  <INSTR_SHP>  <INLU_GRID>  <INMUN_SHP> 
<OUT_GRID> 

&if  [NULL  %outgrid%]  =  .TRUE.  &then 

Sreturn  Sinform  Usage:  ASTRGRD  <INSTR_SHP>  <INLU_GRID>  <INMUN_SHP> 
<OUT_GRID> 

&if  A  [exists  %incover% . shp  -file]  Sthen 

Sreturn  ERROR!  The  SHAPEFILE  %incover%  does  not  exist. 

Sif  A  [exists  %lugrid%  -grid]  Sthen 

Sreturn  ERROR!  The  GRID  %lugrid%  does  not  exist. 

Sif  A  [exists  %inmun%.shp  -file]  Sthen 

Sreturn  ERROR!  The  SHAPEFILE  %inmun%  does  not  exist. 

Sif  [exists  %incover%  -COVER]  or  [exists  %incover%  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  %incover%  already  exist. 

Sif  [exists  astrgrd  -COVER]  or  [exists  astrgrd  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  astrgrd  already  exist. 

Sif  [exists  astrgrdl  -COVER]  or  [exists  astrgrdl  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  astrgrdl  already  exist. 

Sif  [exists  astrgrd2  -COVER]  or  [exists  astrgrd2  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  astrgrd2  already  exist. 

Sif  [exists  astrgrd3  -COVER]  or  [exists  astrgrd3  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  astrgrd3  already  exist. 

Sif  [exists  astrgrd4  -COVER]  or  [exists  astrgrd4  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  astrgrd4  already  exist. 

Sif  [exists  %outgrid%  -COVER]  or  [exists  %outgrid%  -GRID]  Sthen 
Sreturn  ERROR!  The  COVERAGE  or  GRID  %outgrid%  already  exist. 

I  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  converting  road  shapefile  to  cover  * 
/************************************** 

SHAPEARC  %incover% . shp  astrgrdl  order 
ROUTEARC  astrgrdl  order  astrgrd2 
KILL  astrgrdl 

/***************************** 

/*  converting  the  projection  * 

/  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

PROJECT  COVER  astrgrd2  %incover% 
input 

projection  geographic 

units  dd 

datum  nad83 

spheroid  grs80 

parameters 

output 

projection  alters 
units  meters 
datum  nad83 
spheroid  grs80 
parameters 
29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 
0 

end 


KILL  astrgrd2  ALL 
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/****************************************** 

/*  convert  municipal  boundary  to  coverage  * 
/****************************************** 

SHAPEARC  %inmun%.shp  astrgrd4 
CLEAN  astrgrd4 

PROJECT  COVER  astrgrd4  astrgrd5 
input 

projection  geographic 
units  dd 
datum  nad83 
spheroid  grs80 
parameters 

output 

projection  alters 
units  meters 
datum  nad83 
spheroid  grs80 
parameters 

29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 

0 

end 

KILL  astrgrd4 
CLEAN  astrgrd5 

/********************************************** 

/*  reselect  the  A6,  A1  -  A3  from  the  coverage  * 

j  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

RESELECT  %incover%  astrgrdl  LINE 

res  cfcc2  =  ' Alf  OR  cfcc2  =  '  A2'  OR  cfcc2  =  'A3'  OR  cfcc  =  'A63' 

n 

n 

/★★•A-********************************** 

/*  convert  the  road  coverage  to  grid  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

LINEGRID  astrgrdl  astrgrd2 

30 

y 

NODATA 

GRID 

astrgrd  =  (astrgrd2  *  0)  +24 

Q 

KILL  astrgrd2 

/************************************************* 

/*  create  70m  GIRD  along  the  new  street  coverage  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

BUFFER  astrgrdl  astrgrd2  #  #  70  #  LINE 
KILL  astrgrdl  ALL 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  clip  the  buffer  with  municipal  boundary  * 
/******************************************* 

ERASE  astrgrd2  astrgrd5  astrgrdl 

KILL  astrgrd2 

CLEAN  astrgrdl  astrgrd2 
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KILL  astrgrdl 
KILL  astrgrd5 

I  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  convert  to  grid  * 

/******************* 

POLYGRID  astrgrd2  astrgrdl  inside 

30 

Y 

KILL  astrgrd2  ALL 
GRID 

astrgrd2  =  select (astrgrdl ,  'value  =  100') 

KILL  astrgrdl  ALL 

/********************************** 

/*  separate  the  NLCD  landuse  grid  * 

/  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

astrgrd3  =  select (%lugrid%,  'value  =  21  OR  value  =22') 
astrgrd4  =  select (%lugrid%,  'value  =23') 
astrgrdl  =  merge (astrgrd3,  astrgrd2,  astrgrd4) 

KILL  astrgrd2  ALL 
KILL  astrgrd3  ALL 
KILL  astrgrd4  ALL 

astrgrd2  =  reclass (astrgrdl ,  astrgrdrp.txt,  NODATA) 
astrgrd3  =  reclass ( %lugrid% ,  astrgrdrp2.txt,  DATA) 
astrgrd4  =  merge (astrgrd2 ,  astrgrd3) 

KILL  astrgrdl  ALL 
KILL  astrgrd2  ALL 
KILL  astrgrd3  ALL 

%outgrid%  =  merge (astrgrd,  astrgrd4) 

KILL  astrgrd4  ALL 
KILL  astrgrd  ALL 
Q 

KILL  %incover% 

3munbuf . ami 

/********************************************************************* 

/*  3munbuf.aml  * 

/  -k  -k 

/*  This  AML  will  import  the  shape  file  and  * 

/*  change  it  to  Cover/GRID  with  appropriate  projection  * 

/*  then  create  1.5mil  buffer  with  it  * 

/*  * 

/*  Made  by  Yong  Wook  Kim  on  Jan  7,  2003  * 

/*  Modified  by  Woonsup  Choi  on  June  25,  2003  * 

/  ' k'k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k 
j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  setup  the  input  parameter  * 
y***************************** 

&sv  input_exchangein  [response  'Enter  a  municipal  boundary  shape  file  name 
without  extension ' ] 

&sv  input_exchange  =  [locase  %input_exchangein% ] 

&sv  input_bound  [response  'Enter  the  boundary  GRID  name'] 

&sv  input_bnd_exchange  [response  'Enter  the  BOUNDARY  cover  name'] 

&sv  x  [response  'Enter  the  left  end  coordinate  of  the  boundary  grid'] 

&sv  y  [response  'Enter  the  bottom  end  coordinate  of  the  boundary  grid'] 

&sv  nrows  [response  'Enter  the  #  of  rows  of  the  boundary  grid'] 
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Ssv  ncols  [response  'Enter  the  #  of  columns  of  the  boundary  grid'] 

/*  &if  %input_type%  ne  c  and  %input_type%  ne  g  &then 

/*  Sreturn  ERROR!  You  have  to  type  either  ' c'  for  Coverage  or  ' g'  for  GRID 
for  the  type. 

/*  &if  %input_type%  =  c  &then 
&sv  input_type  =  COVER 
/*  &if  %input_type%  =  g  &then 
/*  &sv  input_type  =  GRID 

/*  &sv  change_grid  =  n 

/*  &if  %input_type%  =  COVER  &then 
/*  &do 

/*  &sv  change_grid  [response  'Do  you  want  to  change  the  output  buffer  to 
GRID?  (y/n)'] 

/*  &sv  change_grid  =  [locase  %change_grid% ] 

/*  &if  %change_grid%  ne  y  and  %change_grid%  ne  n  &then 

/*  Sreturn  ERROR!  You  have  to  type  either  'y'  or  'n' 

/*  Send 

&TYPE 

SType  All  THE  EXISTING  COVERAGE  OR  GRID  NAMED  " %input_exchange% , " 
"%input_exchange%l, "  "%input_exchange%buf , "  "%input_bound%l , " 
"%input_bound%2 "  WILL  BE  DELETED  ! ! ! 

&TYPE 

&sv  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

&sv  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  n  and  %delete_yesno%  ne  y  Sthen 
Sreturn  ERROR!  You  have  to  type  either  'y'  or  'n' 

Sif  %delete_yesno%  =  n  sthen 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

/*  Ssv  pro j ect_yesnoin  [response  'Do  you  want  to  perform  the  projection 
change  (y/n) ?  '  ] 

/*  Ssv  proj ect_yesno  =  [locase  %project_yesnoin%] 

/*  Ssv  input_unitin  [response  'Is  the  unit  of  input  file  feet  or  meter? 
(f/m)'] 

/*  Ssv  input_unit  =  [locase  %input  unitin%] 

/*  Sif  %input_unit%  ne  f  and  %input_unit%  ne  m  Sthen 

/*  Sreturn  ERROR!  You  have  to  type  either  'f'  of  'm'  for  the  unit. 

Sif  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

kill  %input_exchange%  all 

Sif  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

Sif  [exists  %input_exchange%buf  -COVER]  or  [exists  %input_exchange%buf  - 
GRID]  Sthen 

kill  %input_exchange%buf  all 

/*Sif  [exists  %input_bound%  -COVER]  or  [exists  %input_bound%  -GRID]  Sthen 
/*  kill  %input_bound%  all 

Sif  [exists  %input_bound%l  -COVER]  or  [exists  %input_bound%l  -GRID]  Sthen 
kill  %input_bound%l  all 

Sif  [exists  %input_bound%2  -COVER]  or  [exists  %input_bound%2  -GRID]  Sthen 
kill  %input_bound%2  all 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  Import  boundary  exchange  file  as  a  cover  and  covert  to  a  grid  * 
/****************************************************************** 
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/*import  cover  %input_bnd_exchange% . eOO  %input_bound% 

/*polygrid  %input_bound%  %input_bound%l 

/*30 

/*y 

/‘rename  %input_bound%l  %input_bound% 

/******************************************** 

/*  changing  a  shape  file  to  poly  Coverage  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

shapearc  %input  exchange% . shp  %input  exchange%  subclass 
clean  %input  exchange%  %input  exchange%l 

kill  %input  exchange% 

rename  %input  exchange%l  %input  exchange% 

/****************************** 

/*  projecting  the  output  file  * 
/****************************** 

pro j ectdef ine  %input_type%  %input_exchange% 
projection  geographic 
datum  nad83 

spheroid  grs80 
units  dd 
parameters 

project  %input  type%  %input  exchange%  %input  exchange%l 
output 

projection  albers 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 
29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 
0 

end 

kill  %input_exchange%  all 

rename  %input_exchange%l  %input  exchange% 
build  %input_exchange% 

/*  Send 

/**************************** 

/*  create  a  1.5  mile  buffer  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  &if  %input_type%  =  COVER  &then 
/*  &do 

buffer  %input_exchange%  %input_exchange%buf  #  #  2414.016 
/*  Send 

/*  &if  %input_type%  =  GRID  &then 
/*  &do 

/*  Send 

/******************************* 

/*  change  the  coverage  to  GRID  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  &if  %change_grid%  =  y  &then 
/*  &do 
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/*  kill  %input_exchange%  all 

rename  %input  exchange%buf  %input  exchange%l 
polygrid  %input  exchange%l  %input  exchange%buf 
30 
n 

-oX-5  ,  'oY'o 

%nrows%,  %ncols% 

kill  %input_exchange%l  all 
/*  Send 


j-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  create  temporary  input  boundary  GRID  * 

I  '■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

COPY  %input_bound%  %input_bound%l 
GRID 

%input_bound%2  =  reclass ( %input_bound%l , boundremap2 ) 
kill  %input_bound%l  all 

J  **•*••*•**•*•***•*••*••*••*••*••*•**•*••*• 

/*  reclassify  the  GRID  * 

j-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%input_exchange%l  =  reclass (%input_exchange%buf ,munbndremap) 
kill  %input  exchange%buf 

/'k-k-k-k-k-k-k'k'k-k-k-k-k'k-k 

/*  merge  GRIDS  * 

/  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%input  bound%l  =  merge (%input  exchange%l , %input  bound%2) 
kill  %input  bound%2 
kill  %input  exchange%l 

/***************** 

/*  clipping  GRID  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

gridclip  %input_bound%l  %input_bound%2  cover  %input_bnd_exchange% 
kill  %input_bound%l  all 

q 

rename  %input  bound%2  %input  exchange%buf 

gridascii  %input  exchange%buf  %input  exchange%buf . asc 

4travel . ami 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  4travel.aml  * 

/*  * 

/*  This  AML  will  produce  Travelminute_Map  GRID  by  using  * 

/*  Landuse  GRID  and  roadmap  exchange  file  * 

/*  Projection  will  be  changed  too,  upon  request  * 

J  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  10,  2003  * 

/*  Modified  by  Woonsup  Choi  on  Feb  19,  2003  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/***************************** 

/*  setup  the  input  parameter  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&sv  input_exchange  [response  'Enter  the  ROAD  cover  exchange  file  w/o 
extension ' ] 

&sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  of  the  road  exchange 
file'] 
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&sv  input_lu  [response  'Enter  the  name  of  the  LANDUSE  GRID'] 

&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  GRID'] 

/*  &sv  buf f er_yesnoin  [response  'Do  you  want  to  buffer  the  limited-access? 
(y/n)  '  ] 

/*  &sv  buffer_yesno  =  [locase  %buf f er_yesnoin% ] 

/*  &if  %buf fer_yesno%  =  y  &then 
/*  &do 

&sv  buffer_item  [response  'Enter  the  class  of  the  limited- 
access  (buffer) ' ] 

& sv  ramp_item  [response  'Enter  the  class  of  the  ramp'] 

/*  Send 


&TYPE 

SType  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 

"%input  exchange%,"  "%input  exchange%l,"  "%input  exchange%2," 
"%input_exchange%3, "  "%input_exchange%4 , "  "%input_exchange%5, " 
"%input_exchange%9,  "  "tspeed"  WILL  BE  DELETED  ! ! ! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

/*  &sv  pro j ect_yesnoin  [response  'Do  you  want  to  perform  the  projection 
change  (y/n) ?  '  ] 

/*  &sv  proj ect^yesno  =  [locase  %project_yesnoin%] 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

/*  Pay  attention  HERE!  This  command  can  be  hidden. 

/*&if  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

/*  kill  %input_exchange%  all 
/* 


&if  [exists  %input  exchange%l  -COVER] 
Sthen 

kill  %input  exchange%l  all 

or 

[exists 

%input  exchange%l 

-GRID] 

&if  [exists  %input  exchange%2  -COVER] 
Sthen 

kill  %input  exchange%2  all 

or 

[exists 

%input  exchange%2 

-GRID] 

&if  [exists  %input  exchange%3  -COVER] 
&then 

kill  %input  exchange%3  all 

or 

[exists 

%input  exchange%3 

-GRID] 

&if  [exists  %input  exchange%4  -COVER] 
&then 

kill  %input  exchange%4  all 

or 

[exists 

%input  exchange%4 

-GRID] 

&if  [exists  %input  exchange%5  -COVER] 
&then 

kill  %input  exchange%5  all 

or 

[exists 

%input  exchange%5 

-GRID] 

&if  [exists  %input  exchange%9  -COVER] 
&then 

or 

[exists 

%input  exchange%9 

-GRID] 

kill  %input_exchange% 9  all 

&if  [exists  tspeed  -COVER]  or  [exists  tspeed  -GRID]  Sthen 
kill  tspeed  all 
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/************************************************************************** 
/*  convert  shapefile  to  coverage  * 

/*  Pay  attention  here!  This  command  can  be  hidden  when  you  use  a  cover,  not 
an  exchange  file  * 

^'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

•k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

/*IMPORT  cover  %input_exchange% . eOO  %input_exchange% 

/*shapearc  %input_exchange% . shp  %input_exchange%  subclass 
/*routearc  %input  exchange%  subclass  %input  exchange%l 
/*kill  %input_exchange% 

/‘rename  %input_exchange%l  %input_exchange% 

y******************************** 

/*  projecting  the  road  coverage  * 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  input_type  =  COVER 

/*  &if  %project_yesno%  =  y  &then 

/*  &do 

I  : k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  give  the  projection  to  the  converted  Coverage  * 
/************************************************* 
pro j ectdef ine  %input_type%  %input_exchange% 
projection  geographic 
datum  nad83 
spheroid  grs80 
units  dd 
parameters 

/***************** 

/*  projecting  it  * 

j  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

project  %input  type%  %input  exchange%  %input  exchange%l 
output 

projection  albers 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 
29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 
0 

end 

kill  %input  exchange% 

RESELECT  %input_exchange% 1  %input_exchange%  line 
res  %value_item%  <=  6  and  %value_item%  >=  1 

n 

n 

kill  %input_exchange%l  all 

y********************************************************** 

/*  buffering  limited-access  hwys  and  converting  to  a  grid  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  &if  %buf fer_yesno%  =  y  &then 
/*  &do 

RESELECT  %input_exchange%  %input_exchange%l  line 
res  %value  item%  =  %buffer  item% 


n 
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n 

BUFFER  %input_exchange%l  %input_exchange%2  #  #  70  0.5  line  flat 
POLYGRID  %input_exchange%2  %input_exchange%4 
30 

Y 

kill  %input_exchange%2  all 

LATTICEPOLY  %input_bound%  %input_exchange%2  ASPECT 
LATTICECLIP  %input_exchange%4  %input_exchange%2  %input_exchange%3 
kill  %input_exchange%4  all 
kill  %input_exchange%2  all 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  converting  the  interstate  to  GRID  * 
/************************************* 

LINEGRID  %input_exchange%l  %input_exchange%4  %value_item% 

30 

Y 

NODATA 

kill  %input_exchange%l  all 

I  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  converting  the  ramp  to  GRID  * 

/******************************* 

RESELECT  %input_exchange%  %input_exchange%l  line 
res  %value  item%  =  %ramp  item% 

n 

n 

LINEGRID  %input_exchange%l  %input_exchange%5  %value_item% 

30 

Y 

NODATA 

kill  %input_exchange%l  all 

I  • k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  converting  the  roads  to  grid  * 

/************************************************ 

LINEGRID  %input_exchange%  %input_exchange%l  %value_item% 

30 

Y 

NODATA 

/*  kill  %input_exchange%  all 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  merge  road  GRID  and  interstate  buffer  * 
/***************************************** 

GRID 

/*  buffer 

%input_exchange%2  =  reclass (%input_exchange%3 , travelremapl ) 
kill  %input_exchange%3  all 

/*  roads 

%input_exchange%9  =  int (%input_exchange%l) 
kill  %input_exchange%l  all 

%input_exchange%l  =  reclass (%input_exchange%9 , travelremap2 ) 
kill  %input_exchange% 9  all 

/*  limited-access  hwy 

%input_exchange%3  =  reclass (%input_exchange%4 , travelremap3) 
kill  %input_exchange%4  all 

/*  ramp 

%input_exchange%4  =  reclass (%input_exchange%5 , travelremap4 ) 
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kill  %input  exchange%5  all 
%input  exchange%5  = 

merge  (%input  exchange%3 , %input  exchange%4 , %input  exchange%2 , %input  exchange% 
'■) 

kill  %input_exchange%l  all 
kill  %input_exchange%2  all 
kill  %input_exchange%3  all 
kill  %input_exchange%4  all 

/*  landuse 

%input_exchange%l  =  reclass (%input_lu%, travelremap5 ) 

%input  exchange%2  =  merge (%input_exchange%5, %input  exchange%l) 
kill  %input_exchange%5  all 
kill  %input_exchange%l  all 

%input_exchange%l  =  1000  /  (26.8224  *  %input_exchange%2 ) 
rename  %input_exchange%2  tspeed 
rename  %input  exchange%l  %output_name% 

q 

/*  %output_name%  =  int (%input_exchange% ) 

/*  kill  %input_exchange% 


6roadatr . ami 

/  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  6roadatr.aml  * 

J  ~k  -k 

/*  This  AML  will  produce  CT  ROAD  ATTRACTOR  MAP  GRID  * 

/*  by  using  roadmap  exchange  file  and  Travel  Minute  Map  * 

/*  Projection  will  be  changed  too  * 

/*  * 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

/*  Modified  by  Woon-sup  Choi  on  Jul  11,  2003  * 

I'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/***************************** 

/*  setup  the  input  parameter  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  x  [response  'Enter  the  left  end  coordinate  of  the  boundary  grid'] 

&sv  y  [response  'Enter  the  bottom  end  coordinate  of  the  boundary  grid'] 

&sv  nrows  [response  'Enter  the  #  of  rows  of  the  boundary  grid'] 

&sv  ncols  [response  'Enter  the  #  of  columns  of  the  boundary  grid'] 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&sv  input_exchange  [response  'Enter  the  ROAD  EXCHANGE  FILE  name  without 
extension ' ] 

&sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  of  the  road  exchange 
file'] 

&sv  road_class  [response  'Enter  the  CLASS  of  ROADS  and  STREETS'] 

&sv  travel_map  [response  'Enter  the  name  of  the  TRAVEL  MINUTE  MAP'] 

&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  grid'] 

&TYPE 

SType  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%input_exchange% , "  "%input_exchange%l, "  "%input_exchange%2 "  WILL  BE  DELETED 
!  !  ! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 
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Sreturn  ERROR!  You  have  to  type  either  ' y'  of  ' n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

kill  %input_exchange%  all 

&if  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

&if  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 

j  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  convert  shapefile  coverage  * 

/****************************** 

IMPORT  cover  %input_exchange% . eOO  %input_exchange%x 
RESELECT  %input_exchange%x  %input_exchange%l  line 
res  %value_item%  =  %road_class% 

n 

n 

kill  %input  exchange%x  all 

RESELECT  %input_exchange% 1  %input_exchange%  line 
res  fetype  lk  'rd' 

n 

n 

kill  %input_exchange%l  all 

/******************************** 

/*  projecting  the  road  coverage  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  input_type  =  COVER 

/************************************************* 

/*  give  the  projection  to  the  converted  Coverage  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

pro j ectdef ine  %input_type%  %input_exchange% 

projection  geographic 

datum  nad27 

spheroid  clarkel866 

units  DD 

parameters 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  projecting  it  * 

/***************** 

project  %input  type%  %input  exchange%  %input  exchange%l 
output 

projection  alters 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 
29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 
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0 

end 

kill  %input_exchange%  all 

/************************************************ 

/*  converting  the  roads  to  grid  and  clipping  it  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

LINEGRID  %input_exchange%l  %input_exchange%  %value_item% 

30 

n 

"OX'S  f  'oY'o 

%nrows%,  %ncols% 

NODATA 

kill  %input_exchange%l  all 
grid 

/**************************** 

/*  calculating  costdistance  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%input_exchange%l  =  costdistance (%input_exchange%, %travel_map%, #, #, 120, #) 
%output_name%l  =  merge  (int (%input_exchange%l ) , %input_bound%  *  0  +  120) 
kill  %input_exchange%l  all 
kill  %input_exchange%  all 

q 

rename  %output  name%l  %output_name% 
gridascii  %output_name%  %output_name% . asc 

9nogrow . ami 

y*************************************************************** 

/*  9nogrow.aml  * 

J  -k  -k 

/*  This  AML  will  produce  NO  GROWTH  ZONE  MAP  GRID  * 

/*  by  using  various  exchange  files  such  as  natrual  areas,  * 

/*  nature  preserves,  state  conservation  areas,  and  so  on  * 

/*  Projection  will  be  changed  too  * 

/  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

/*  Modified  by  Woonsup  Choi  Jul  11,  2003  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

y***************************** 

/*  setup  the  input  parameter  * 
y***************************** 

&sv  x  [response  'Enter  the  left  end  coordinate  of  the  boundary  grid'] 

&sv  y  [response  'Enter  the  bottom  end  coordinate  of  the  boundary  grid'] 

&sv  nrows  [response  'Enter  the  #  of  rows  of  the  boundary  grid'] 

&sv  ncols  [response  'Enter  the  #  of  columns  of  the  boundary  grid'] 

&sv  pub_exchange  [response  'Enter  the  PUBLIC  LAND  SHAPE  FILE  name  without 
extension ' ] 

/*  Before  running,  combined  all  the  protected  area  shapefiles  into  one  file 
in  geographic  coord. 

/*&sv  np_exchange  [response  'Enter  the  NATURAL  PRESERVES  EXCHANGE  FILE  name 
without  extension  '  ] 

/*&sv  sc^exchange  [response  'Enter  the  STATE  CONSERVATION  AREAS  EXCHANGE 
FILE  name  without  extension'] 

/*&sv  sf_exchange  [response  'Enter  the  STATE  FORESTS  EXCHANGE  FILE  name 
without  extension'] 

/*&sv  sfw_exchange  [response  'Enter  the  STATE  FISH  EXCHANGE  FILE  name 
without  extension'] 

/*&sv  sp_exchange  [response  'Enter  the  STATE  PARK  EXCHANGE  FILE  name  without 
extension  '  ] 

&sv  rd_exchange  [response  'Enter  the  ROAD  EXCHANGE  FILE  name  without 
extension  '  ] 
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&sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  of  the  road  exchange 
file'] 

&sv  road_class  [response  'Enter  the  CLASS  of  the  LIMITED-ACCESS  HWY ' ] 

&sv  input_grid  [response  'Enter  the  name  of  the  BOUNDARY  grid'] 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&TYPE 

&Type  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%pub_exchange% "  "%rd_exchange%"  "%rd_exchange%l"  WILL  BE  DELETED  ! ! ! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %pub_exchange%  -COVER]  or  [exists  %pub_exchange%  -GRID]  Sthen 
kill  %pub_exchange%  all 

&if  [exists  %pub_exchange%l  -COVER]  or  [exists  %pub_exchange%l  -GRID]  Sthen 
kill  %pub_exchange%l  all 

/*&if  [exists  %np_exchange%  -COVER]  or  [exists  %np_exchange%  -GRID]  Sthen 
/*  kill  %np_exchange%  all 

/*&if  [exists  %sc_exchange%  -COVER]  or  [exists  %sc_exchange%  -GRID]  Sthen 
/*  kill  %sc_exchange%  all 

/*&if  [exists  %sf_exchange%  -COVER]  or  [exists  %sf_exchange%  -GRID]  &then 
/*  kill  %sf_exchange%  all 

/*&if  [exists  %sfw_exchange%  -COVER]  or  [exists  %sfw_exchange%  -GRID]  &then 
/*  kill  %sfw_exchange%  all 

/*&if  [exists  %sp_exchange%  -COVER]  or  [exists  %sp_exchange%  -GRID]  &then 
/*  kill  %sp_exchange%  all 

&if  [exists  %rd_exchange%  -COVER]  or  [exists  %rd_exchange%  -GRID]  &then 
kill  %rd_exchange%  all 

&if  [exists  %rd_exchange%l  -COVER]  or  [exists  %rd_exchange%l  -GRID]  &then 
kill  %rd_exchange%l  all 

&if  [exists  %rd_exchange%2  -COVER]  or  [exists  %rd_exchange%2  -GRID]  &then 
kill  %rd_exchange%2  all 

&if  [exists  %input_grid%l  -COVER]  or  [exists  %input_grid%l  -GRID]  &then 
kill  %input  grid%l  all 

/******** 

/*  road  * 

/•k-k-k-k-k-k-k-k 

/*  import  coverage 

IMPORT  cover  %rd_exchange% . eOO  %rd_exchange%x 

/*  extracting  the  interstate  * 

RESELECT  %rd_exchange%x  %rd_exchange%  line 
res  %value_item%  =  %road_class% 

n 

n 

kill  %rd_exchange%x  all 

pro j ectdef ine  cover  %rd_exchange% 

projection  geographic 

datum  nad27 

spheroid  clarkel866 

units  DD 

parameters 
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/*  projecting  it 

project  cover  %rd_exchange%  %rd_exchange%l 
output 

projection  alters 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 

29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 

0 

end 

kill  %rd_exchange%  all 

/*  create  70  meter  buffer  to  interstate  * 

BUFFER  %rd_exchange%l  %rd_exchange%  #  #  70  0.002 
POLYGRID  %rd_exchange%  %rd_exchange%2 

30 
n 

SXS  f  'oY'o 

%nrows%,  %ncols% 

kill  %rd_exchange% 
kill  %rd_exchange%l 

rename  %rd_exchange%2  %rd_exchange% 

y*************************************************************************** 

* 

/*  convert  shapefile  to  coverage  with  proper  projection  and  convert  to  GRID 

•k 

■k 

&sv  input_type  =  COVER 

/************************ 

/*  Any  public  land  area  * 
y************************ 

SHAPEARC  %pub_exchange% . shp  %pub_exchange%  subclass 
clean  %pub_exchange%  %pub_exchange%l 
kill  %pub_exchange% 

rename  %pub_exchange%l  %pub_exchange% 

/*  project  define 

pro j ectdef ine  %input_type%  %pub_exchange% 
projection  geographic 
datum  nad83 
spheroid  grs80 
units  dd 
parameters 

/*  projecting  it 

project  %input_type%  %pub_exchange%  %pub_exchange%l 
output 

projection  albers 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 
29  30  00 
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45  30  00 
-96  00  00 
23  00  00 
0 
0 

end 

kill  %pub_exchange%  all 
build  %pub_exchange%l 

/*  convert  to  grid 

POLYGRID  %pub_exchange%l  %pub_exchange% 

30 

n 

-OX'S  r  -sy  -o 

%nrows%,  %ncols% 

kill  %pub_exchange%l  all 

grid 

%pub_exchange%l  =  %pub_exchange%  *  0  t  1 
kill  %pub_exchange% 

q 

rename  %pub_exchange%l  %pub_exchange% 

/********************** 

/*  merging  final  GRID  * 

/********************** 

grid 

%rd_exchange%l  =  %rd_exchange%  *  0  t  1 
%rd_exchange%2  =  %input_grid%  *  0 

%output_name%l  =  merge (%rd_exchange%l , %pub_exchange%, %rd_exchange%2 ) 
%output  name%  =  %output  name%l  *  %input  grid% 

q 

gridascii  %output_name%  %output_name% . asc 

/*kill  %na_exchange%  all 

/*kill  %np_exchange%  all 

/*kill  %sc_exchange%  all 

/*kill  %sf_exchange%  all 

/*kill  %sfw_exchange%  all 

kill  %rd_exchange%  all 

kill  %rd_exchange%l  all 

kill  %rd_exchange%2  all 

kill  %pub_exchange%  all 

kill  %output_name%l  all 

lOrampatr . ami 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  lOrampatr . ami  * 

/  *  * 

/*  This  AML  will  produce  RAMP  ATTRACTOR  MAP  GRID  * 

/*  by  using  roadmap  exchange  file  and  Travel  Minute  Map  * 

/*  Projection  will  be  changed  too  * 

/  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/***************************** 

/*  setup  the  input  parameter  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  x  [response  'Enter  the  left  end  coordinate  of  the  boundary  grid'] 
&sv  y  [response  'Enter  the  bottom  end  coordinate  of  the  boundary  grid'] 
&sv  nrows  [response  'Enter  the  #  of  rows  of  the  boundary  grid'] 

&sv  ncols  [response  'Enter  the  #  of  columns  of  the  boundary  grid'] 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 
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Ssv  input_exchange  [response  'Enter  the  ROAD  EXCHANGE  FILE  name  without 
extension ' ] 

&sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  of  the  road  exchange 
file'] 


Ssv 

road 

class 

[response 

' Enter 

Ssv 

travel  map 

[response 

' Enter 

Ssv 

input 

grid 

[response 

' Enter 

the  CLASS  of  the  RAMP'] 

the  name  of  the  TRAVEL  MINUTE  MAP'] 

the  name  of  the  BOUNDARY  grid'] 


&TYPE 

SType  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%input_exchange% ,  "  "%input_exchange%l , "  "%input_exchange%2 "  WILL  BE  DELETED 

i  i  i 

&TYPE 


Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n)  '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

kill  %input_exchange%  all 

&if  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

&if  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 
Sif  [exists  xxx  -COVER]  or  [exists  xxx  -GRID]  Sthen 
kill  xxx  all 


I : kkkkkkkkkkkkkkkkkkk 

/*  import  coverage  * 

/******************* 

IMPORT  cover  %input_exchange% . eOO  %input_exchange%x 
RESELECT  %input_exchange%x  %input_exchange%l  line 
res  %value_item%  =  %road_class% 

n 

n 

kill  %input  exchange%x  all 

/******************************** 

/*  projecting  the  ramp  coverage  * 

^-k-k-k-k-k-k-kie-k-k-k-k-k-k-k-k-k-k-kie-k-k-k-k-k-k-k-k-k-k-k-k 

Ssv  input_type  =  COVER 

/kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

/*  give  the  projection  to  the  converted  Coverage  * 

/kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

pro j ectdef ine  %input  type%  %input  exchange%l 

projection  geographic 

datum  nad27 

spheroid  clarkel866 

units  DD 
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parameters 


I : k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

/*  projecting  it  * 

/***************** 

project  %input  type%  %input  exchange%l  %input_exchange% 
output 

projection  alters 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 

29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 

0 

end 

kill  %input_exchange%l  all 

/************************************************ 

/*  converting  the  roads  to  grid  and  clipping  it  * 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

LINEGRID  %input_exchange%  %input_exchange%l  %value_item% 

30 
n 

'oX'o  f  'oY'o 

%nrows%,  %ncols% 

NODATA 

kill  %input_exchange%  all 
grid 

rename  %input  exchange%l  %input  exchange% 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  calculating  costdistance  * 

/**************************** 

%input_exchange%l  =  costdistance (%input_exchange%, %travel_map%, #, #, 120, #) 
%output_name%  =  merge  ( int ( %input_exchange%l ) , %input_gr id%  *  0  +  120) 
kill  %input_exchange%  all 
kill  %input_exchange%l  all 

q 

gridascii  %output_name%  %output_name% . asc 
llrdintatr. ami 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  llrdintatr . ami  * 

/  *  -k 

/*  This  AML  will  produce  RD  INTERSETION  MAP  GRID  * 

/*  by  using  roadmap  exchange  file  and  Travel  Minute  Map  * 

/*  Projection  will  be  changed  too  * 

/  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

/*  Modified  by  Woonsup  Choi  on  July  15,  2003  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/***************************** 

/*  setup  the  input  parameter  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 
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/*&sv  input_exchange  [response  'Enter  the  INTERSECTION  POINT  EXCHANGE  FILE 
name  without  extension'] 

&sv  input_exchange  [response  'Enter  the  INTERSECTION  POINT  SHAPE  FILE  name 
without  extension  '  ] 

/*  A  node  shapefile  must  be  created  before  running  this  AML.  Do  it  with  the 
" intersection_point . apr "  project  file.  Select  only  Node_type  =  "real" 

&sv  travel_map  [response  'Enter  the  name  of  the  TRAVEL  MINUTE  MAP'] 

&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  COVER'] 

&sv  input_grid  [response  'Enter  the  name  of  the  BOUNDARY  grid'] 

&sv  x  [response  'Enter  the  left  end  coordinate  of  the  boundary  grid'] 

&sv  y  [response  'Enter  the  bottom  end  coordinate  of  the  boundary  grid'] 

&sv  nrows  [response  'Enter  the  #  of  rows  of  the  boundary  grid'] 

&sv  ncols  [response  'Enter  the  #  of  columns  of  the  boundary  grid'] 

&TYPE 

SType  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%input_exchange% , "  "%input_exchange%, "  "%input_exchange%l, " 
"%input~exchange%2"  WILL  BE_ DELETED  !!! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 
Sthen 

kill  %input_exchange%  all 

&if  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

&if  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 

/******************* 

/*  import  coverage  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*IMPORT  cover  %input_exchange% . eOO  %input_exchange% 

/******************************************* 

/*  Converting  the  shapefile  to  a  cover  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

shapearc  %input  exchange% . shp  %input  exchange% 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  Projecting  the  cover  * 

/****************************************** 
pro j ectdef ine  cover  %input_exchange% 
projection  geographic 
datum  nad27 
spheroid  clarkel866 
units  DD 
parameters 

project  cover  %input  exchange%  %input  exchange%l 
output 
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projection  alters 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 

29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 

0 

end 

kill  %input  exchange% 

rename  %input  exchange%l  %input  exchange% 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 
/*  Changing  the  point  cover  to  grid  * 
/************************************* 

pointgrid  %input  exchange%  %input  exchange%l 

30 
n 

"OX'S  f  'oY'o 

%nrows%,  %ncols% 

NODATA 

kill  %input_exchange%  all 
GRID 

/*GRIDCLIP  %input_exchange%l  %input_exchange%  cover  %input_bound% 

/*kill  %input_exchange%l  all 

rename  %input  exchange%l  %input  exchange% 

/**************************** 

/*  calculating  costdistance  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%input_exchange%l  =  costdistance (%input_exchange%, %travel_map%, #, #, 120, #) 
%output_name%l  =  merge ( int ( %input_exchange%l ), %input_grid%  *  0  +  120) 
gridclip  %output  name%l  %output  name%  cover  %input  bound% 

kill  %input_exchange%  all 
kill  %input_exchange%l  all 
kill  %output_name%l  all 

q 

gridascii  %output_name%  %output_name% . asc 

/*  You  must  have  merged  all  the  necessary  DEM  quadrangles  that  cover  the 
study  area.  It's  in  geographic  coord. 

/*  You  must  have  a  projection  file  "geog_alb.txt"  in  the  workspace. 
/***************************** 

/*  setup  the  input  parameter  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  clip_out  [response  'Type  the  name  of  the  output  SLOPE  GRID  name'] 

&sv  demgrid  [response  'Enter  a  merged  dem  name'] 

&sv  clip_bound  [response  'Type  the  name  of  the  boundary  Cover'] 

/*&sv  index  =  0 

/*&sv  biglat  =  %demprefix%0 


&TYPE 
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SType  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%clip_out%, "  " %clip_out%l , 
WILL  BE  DELETED  ! ! ! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  okay  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 


Sreturn 

ERROR ! 

You  have  to  type 

either 

'  y '  of 

'n'  for  the  answer. 

&if 

%delete 

yesno! 

b  =  n 

Sthen 

Sreturn 

Please 

take 

care  of 

the 

file  name  and 

relaunch  the  script 

&if 

[exists 

%clip 

out% 

-COVER] 

or  [ 

exists 

%clip  out%  -GRID]  Sthen 

kill  %clip  out 

%  all 

&if 

[exists 

%clip 

out%l 

-COVER] 

or 

[exists 

%clip 

out%l  -GRID]  Sthen 

kill  %clip_out%l  all 

/****************************** 

/*  projecting  the  merged  grid  * 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

project  grid  %demgrid%  %demgrid%l  geog_alb.txt  bilinear  30 
kill  %demgrid% 


grid 

I  • k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  clip  the  image  GRID  by  using  the  boundary  cover  * 
/*************************************************** 
&if  [exists  %clip_bound%  -cover]  &then 
&do 

gridclip  %demgrid%l  %demgrid%  cover  %clip  bound% 
Send 

kill  %demgrid%l 

j  ' k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  create  slop_map  * 

/******************* 

%clip_out%  =  slope (%demgrid%) 

%clip_out%l  =  int  (%clip_out%) 
kill  %clip_out%  all 
rename  %clip_out%l  %clip  out% 
kill  %clip_out%l  all 

q 

gridascii  %clip_out%  %clip_out% . asc 


13stateatr . ami 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  13stateatr . ami  * 

/*  * 

/*  This  AML  will  produce  ST  RD  ATTRACTOR  MAP  GRID  * 

/*  by  using  roadmap  exchange  file  and  Travel  Minute  Map  * 

/*  Projection  will  be  changed  too  * 

J  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

/*  Modified  by  woonsup  Choi  on  Jul  10,  2003  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 


/*  setup  the  input  parameter  * 

J  •k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


&SV 

&SV 

&SV 

&SV 


x  [response  'Enter  the 
y  [response  'Enter  the 
nrows  [response  'Enter 
ncols  [response  'Enter 


left  end  coordinate  of  the  boundary  grid'] 
bottom  end  coordinate  of  the  boundary  grid'] 
the  #  of  rows  of  the  boundary  grid'] 
the  #  of  columns  of  the  boundary  grid' ] 


&sv  input_exchange  [response  'Enter  the  ROAD  EXCHANGE  FILE  name'] 
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/*  This  must  be  a  cover  converted  from  TIGER  shapefile  in  geographic 
coordinate 

&sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  of  the  road  exchange 
file'] 

&sv  road_class  [response  'Enter  the  CLASS  of  the  STATE  HIGHWAY'] 

&sv  travel_map  [response  'Enter  the  name  of  the  TRAVEL  MINUTE  MAP'] 

&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  grid'] 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&TYPE 

&Type  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%input_exchange% , "  "%input_exchange%l, "  "%input_exchange%2 "  WILL  BE  DELETED 
!  !  ! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %input_exchange%  -COVER]  or  [exists  %input_exchange%  -GRID] 

Sthen 

kill  %input_exchange%  all 

&if  [exists  %input_exchange%l  -COVER]  or  [exists  %input_exchange%l  -GRID] 
Sthen 

kill  %input_exchange%l  all 

&if  [exists  %input_exchange%2  -COVER]  or  [exists  %input_exchange%2  -GRID] 
Sthen 

kill  %input_exchange%2  all 

/******************* 

/*  import  coverage  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

IMPORT  cover  %input_exchange% . eOO  %input_exchange%x 
RESELECT  %input_exchange%x  %input_exchange%  line 
res  %value_item%  =  %road_class% 

n 

n 

kill  %input  exchange%x  all 

/‘rename  %input  exchange%l  %input  exchange% 

/******************************** 

/*  projecting  the  ramp  coverage  * 

j  ■ k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

&sv  input_type  =  COVER 

/************************************************* 

/*  give  the  projection  to  the  converted  Coverage  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

pro j ectdef ine  %input_type%  %input_exchange% 

projection  geographic 

datum  nad27 

spheroid  clarkel866 

units  DD 

parameters 
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/***************** 

/*  projecting  it  * 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

project  %input  type%  %input  exchange%  %input  exchange%l 
output 

projection  alters 
datum  nad83 
spheroid  grs80 
units  meters 
parameters 

29  30  00 
45  30  00 
-96  00  00 
23  00  00 
0 

0 

end 

kill  %input_exchange%  all 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 
/*  converting  the  roads  to  grid  and  clipping  it  * 
/************************************************ 

LINEGRID  %input_exchange%l  %input_exchange% 

30 
n 

"OX'S  f  'oY'o 

%nrows%,  %ncols% 

NODATA 

kill  %input_exchange%l  all 
GRID 

y**************************** 

/*  calculating  costdistance  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%input_exchange%l  =  costdistance (%input_exchange%, %travel_map%, #, #, 120, #) 
%output_name%  =  merge (int (%input_exchange%l ) ,  %input_bound%  *  0  +  120) 

kill  %input_exchange%  all 
kill  %input_exchange%l  all 

q 

gridascii  %output_name%  %output_name% . asc 
14treeatr . ami 

y*************************************************************** 

/*  14treeatr . ami  * 

/*  * 

/*  This  AML  will  produce  TREE  ATTRACTOR  MAP  GRID  * 

/*  by  using  LANDUSE  GRID  * 

/*  Projection  will  not  be  changed  * 

J  -k  -k 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003  * 

y*************************************************************** 

y***************************** 

/*  setup  the  input  parameter  * 
y***************************** 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&sv  input_exchange  [response  'Enter  the  LANDUSE  GRID  name'] 

/*  &sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  name'] 
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&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  cover'] 

&TYPE 

&Type  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%output_name%l , "  "%output_name%2 "  WILL  BE  DELETED  !!! 

&TYPE 

Ssetvar  delete^yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %output_name%l  -COVER]  or  [exists  %output_name%l  -GRID]  Sthen 
kill  %output_name%l  all 

&if  [exists  %output_name%2  -COVER]  or  [exists  %output_name%2  -GRID]  Sthen 
kill  %output_name%2  all 

/***************************************** 

/*  select  forest  cells  from  landuse  GRID  * 

/  ' k**************************************** 

/*  & sv  command  =  ' %value_item%  >=  41  and  %value_item%  <=  43' 

/*  Stype  %value_item% 

/*  Stype  %command% 

GRID 

%output_name%l  =  select (%input_exchange%, 'value  >=  41  and  value  <=  43') 

ft k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  create  the  euclidian  distance  * 

/********************************* 

%output_name%  =  eucdistance (%output_name%l ) 
kill  %output_name%l 

%output_name%l  =  int (%output_name%) 
kill  %output_name% 

rename  %output  name%l  %output_name% 

q 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  clipping  the  output  GRID  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

LATTICECLIP  %output_name%  %input_bound%  %output_name%l 
kill  %output_name%  all 
/*kill  %output_name%l  all 
rename  %output  name%l  %output_name% 
gridascii  %output_name%  %output_name% . asc 

15treeatr . ami 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  15treeatr . ami  * 

/  -k  -k 

/*  This  AML  will  produce  WATER  ATTRACTOR  MAP  GRID 
/*  by  using  LANDUSE  GRID  * 

/*  Projection  will  not  be  changed  * 

/*  * 

/*  Made  by  Yong  Wook  Kim  on  Jan  15,  2003 
/*  Modified  by  Woonsup  Choi  on  Jul  10,  2003 
/***************************************************** 


********** 
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/***************************** 

/*  setup  the  input  parameter  * 

/***************************** 

&sv  output_name  [response  'Enter  the  FINAL  OUTPUT  GRID  name'] 

&sv  input_exchange  [response  'Enter  the  LANDUSE  GRID  name'] 

/*  &sv  value_item  [response  'Enter  the  VALUE  ITEM  FIELD  name'] 

&sv  input_bound  [response  'Enter  the  name  of  the  BOUNDARY  cover'] 

&TYPE 

&Type  All  THE  EXISTING  Coverage  or  GRID  NAMED  "%output_name%, " 
"%output_name%l , "  "%output_name%2 "  WILL  BE  DELETED  !!! 

&TYPE 

Ssetvar  delete_yesnoin  [response  'Is  it  oaky  for  you?  (y/n) '  ] 

Ssetvar  delete_yesno  =  [locase  %delete_yesnoin% ] 

&if  %delete_yesno%  ne  y  and  %delete_yesno%  ne  n  &then 

Sreturn  ERROR!  You  have  to  type  either  'y'  of  'n'  for  the  answer. 

&if  %delete_yesno%  =  n  &then 

Sreturn  Please  take  care  of  the  file  name  and  relaunch  the  script. 

&if  [exists  %output_name%  -COVER]  or  [exists  %output_name%  -GRID]  Sthen 
kill  %output_name%  all 

&if  [exists  %output_name%l  -COVER]  or  [exists  %output_name%l  -GRID]  Sthen 
kill  %output_name%l  all 

&if  [exists  %output_name%2  -COVER]  or  [exists  %output_name%2  -GRID]  Sthen 
kill  %output_name%2  all 

/**************************************** 

/*  select  water  cells  from  landuse  GRID  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

GRID 

%output_name%l  =  select (%input_exchange%, 'value  =  11  or  value  =  91  or  value 
=  92  '  ) 

/********************************* 

/*  create  the  euclidian  distance  * 

I  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

%output_name%  =  int (eucdistance (%output_name%l ) ) 
kill  %output_name%l 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

/*  clipping  the  output  GRID  * 

j  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

GRIDCLIP  %output_name%  %output_name%l  cover  %input_bound% 

q 

kill  %output_name%  all 
/*kill  %output_name%l  all 
rename  %output  name%l  %output_name% 
gridascii  %output_name%  %output_name% . asc 
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Appendix  C:  GRASS  Makefile 

############################################################################ 

# 

#  Last  modified:  Tue  Dec  16  22:05:17  CST  2003 

# 

#  This  file  is  used  by  the  unix  make  program  to  create  the 

#  LEAM  landuse  evolution  input  files  using  grass  commands. 

# 

#  To  use, 

#  Start  GRASS  and  use  a  mapset  that  will  contain  the  results 

#  Create  a  directory  named  cell  in  this  mapset  (if  not  already  there) 

#  Copy  this  Makefile  to  that  directory 

#  Set  the  name  of  the  needed  input  maps  (see  below) 

#  Set  the  size  of  the  cities  (see  below) 

#  Set  the  attractor  weights  (see  below) 

#  From  the  cell  directory,  run  make  on  this  file 

############################################################################ 

# 

#  Input  Maps  and  Map  Categories 

#  Make  sure  the  input  maps  use  the  following  categories  . . . 

#  SMALL_CITY  -  Location  of  "small  cities" 

#  0  No  small  city 

#  1  Small  city 

#  MED_CITY  -  Location  of  "medium  cities" 

#  0  No  small  city 

#  1  Small  city 

#  LARGE_CITY  -  Location  of  "large  cities" 

#  0  No  small  city 

#  1  Small  city 

#  NO_GROWTH  -  Explicit  identification  of  areas  of  no  growth 

#  0  No  restrictions 

#  1  No  growth  areas 

#  DEM 

#  Categories  are  meters  above  sea  level 

#  MUN I C I PAL_BOUNDARY 

#  City  boundaries 

#  STUDY_AREA 

#  0  Outside  study  area 

#  1  Inside  study  area 

#  ROADS 

#  0  No  road 

#  1  Limited  access  highway 

#  2  Federal  highway 

#  3  State  highway 

#  4  County  road 

#  5 

#  6  Ramps 

#  7  Private  road 

#  LANDCOVER  -  NLCD  Landcover  Map 

#  11  Open  water 

#  12  Perennial  Ice/Snow 

#  21  Low  Intensity  Residential 

#  22  High  Intensity  Residential 

#  23  Commercial/Industrial/Transportation 

#  31  Bare  Rock/Sand/Clay 

#  32  Quarries/Strip  Mines/Gravel  Pits 

#  33  Transitional 
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# 

41 

Deciduous  Forest 

# 

42 

Evergreen  Forest 

# 

43 

Mixed  Forest 

# 

51 

Shrubland 

# 

61 

Orchards /Vineyards /Other 

# 

71 

Grasslands /Herbaceous 

# 

81 

Pastures/Hay 

# 

82 

Row  Crops 

# 

83 

Small  Grains 

# 

84 

Fallow 

# 

85 

Urban/Recreational  Grasses 

# 

91 

Woody  Wetlands 

# 

92 

Emergent  Herbaceous  Wetlands 

# 

############################################################################ 

# 

#  Provide  names  of  the  input  maps  (name@mapset) 
SMALL_CITY=small_city@LEAMBase 
MED_CITY=med_city@LEAMBase 
LARGE_CITY=large_city@LEAMBase 
NO_GROWTH=br_nogrowth@ LEAMBa se 
DEM=dem@kbrock 

MUN I C I PAL_BOUNDARY=br_mnbnd@  LEAMBa  s e 
STUDY_AREA=Bragg_county@LEAMBase 
LANDCOVER=landcover @LEAMBase 
ROADS =road@ LEAMBa se 

# 

############################################################################ 

# 

#  Set  the  population  average  for  large,  medium,  and  small  cities.  These 

#  will  be  used  to  create  the  cities_attractor  map,  which  combines  these 

#  populations  driving  time  maps. 

LARGE POP=95 000 

MEDIUMPOP=22000 

SMALLPOP=3500 

# 

############################################################################ 

# 

#  Set  the  multipliers  for  residential  growth  drivers.  These  are  # 

multiplied  times  their  respective  0-1  index  maps.  The  results  are  #  summed 
and  divided  by  the  sum  of  the  multipliers  to  give  an  overall  #  0-1  growth 
attractor  index  map. 

CT_RD_ATTRACTOR_RES_MULT= . 7 

DEM_RES_MULT=1 

RAMP_ATTRACTOR_RES_MULT=0 

ST_RD_ATTRACTOR_RES_MULT=2 

TREE_ATTRACTOR_RES_MULT=. 5 

WATER_ATTRACTOR_RES_MULT= . 2 

CITIES_ATTRACTOR_RES_MULT=5 

############################################################################ 

# 

############################################################################ 

# 

############################################################################ 

# 

#  IDEALLY  NO  EDITING  IS  REQUIRED  BELOW  HERE 

############################################################################ 

# 

############################################################################ 

# 

############################################################################ 

# 
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LEmaps=LElandcover  LEwater_att  LEforest_att  LEroad_att  LEramp_att 
LEcity_buff  LEinter sect_att  LEstaterd_att  LEslope  LEno_growth  LEboundary 
LEcities_att 

IndexMaps=ST_RD_ATTRACTOR_RES  ST_RD_ATTRACTOR_OPENSPACE  \ 
ST_RD_ATTRACTOR_COM_IND  DEM_RES  DEM_OPENSPACE  DEM_COM_IND  \ 
WATER_ATTRACTOR_RES  WATER_ATTRACTOR_OPENSPACE  WATER_ATTRACTOR_COM_IND  \ 

T RE E_AT T RAC T  0R_0 PENSPACE  TREE_ATTRACTOR_COM_IND  T RE E_AT T RAC T 0 R_RE S  \ 
RD_INTERSECTION_RES  RD_INTERSECTION_OPENSPACE  RD_INTERSECTION_COM_IND  \ 

RAM P_AT T RAC T 0 R_RE S  RAMP_ATTRACTOR_OPENSPACE  RAM P_AT T RAC T 0 R_C 0M__ I N D  \ 
CT_RD_ATTRACTOR_RES  CT_RD_ATTRACTOR_OPENSPACE  CT_RD_ATTRACTOR_COM_IND  \ 
CITIES_ATTRACTOR_RES  CITIES_ATTRACTOR_OPENSPACE  CITIES_ATTRACTOR_COM_IND  \ 

OtherMaps=distCost  hwyBuff  road  travelTime30  LElandcover  municipalBoundary  \ 
intersection  smallCity  mediumCity  largeCity  dem  travelSpeed30  travelTime90  \ 
largeCityTime  largeCityTime30  largeCityTime90  mediumCityTime  \ 
mediumCityTime30  mediumCityTime90  smallCityTime  smallCityTime30  \ 
smallCityTime90  landTravelSpeed30  roadTravelSpeed30  travelableLandSpeed  \ 
travelSpeed90  travelTime 

LEfiles=$ (shell  echo  ${LEmaps}  |  sed  -e  's/  /,/g') 

Otherfiles=$ ( shell  echo  ${OtherMaps}  |  sed  -e  ' s/  /,/g') 

Indexfiles=$ ( shell  echo  ${IndexMaps}  |  sed  -e  's/  /,/g') 

all:  ${LEmaps} 

indexmaps :  ${IndexMaps} 

###################################################################### 

#  This  is  the  "final"  step,  combining  residential  attractors 

#  into  an  overall  residential  attractor  map 

ATTRACTOR_RES : CT_RD_ATTRACTOR_RES  DEM_RES  RAM P_AT T RAC T 0 R_RE S  \ 
ST_RD_ATTRACTOR_RES  T RE E_AT T RAC T OR_RE S  WAT E R_AT T RAC T 0 R_RE S  \ 
CITIES_ATTRACTOR_RES  LEno_growth  LElandcover  LEboundary 
r.mapcalc  ATTRACTOR_RES=eval ' (  \ 
developable=if (LEno_growth  | |  \ 

LElandcover==ll  | |  LElandcover==2 1  | |  \ 

LElandcover==22  | |  LElandcover==23  | |  LElandcover==24  | |  \ 
LElandcover==91  | |  LElandcover==92  | |  LElandcover==85,  0,  1),  \ 
developable  *  ( C T_R D_AT T RAC T 0 R_RE S  *  $ { CT_RD_ATTRACTOR_RES_MULT }  +  \ 
DEM_RES  *  $ { DEM_RES_MULT }  +  \ 

RAM  P_AT  T  RAC  T  0  R_RE  S  *  $ { RAMP_ATTRACTOR_RES_MULT }  +  \ 

ST_RD_ATTRACTOR_RES  *  $ { ST_RD_ATTRACTOR_RES_MULT }  +  \ 

T RE E_AT T RAC T 0 R__RE S  *  $  {  TREE_ATTRACTOR_RES_MULT }  +  \ 

WATER_ATTRACTOR_RES  *  $ { WATER_ATTRACTOR_RES_MULT }  +  \ 
CITIES_ATTRACTOR_RES  *  $ { CITIES_ATTRACTOR_RES_MULT } )  /  \ 

( $ { CT_RD_ATTRACTOR_RES_MULT }  +  $ { DEM_RES_MULT }  +  \ 

$ { RAMP_ATTRACTOR_RES_MULT }  +  $ { ST_RD_ATTRACTOR_RES_MULT }  +  \ 

$ { TREE_ATTRACTOR_RES_MULT }  +  $ { WATER_ATTRACTOR_RES_MULT }  +  \ 

$ { CITIES_ATTRACTOR_RES_MULT } ) ) ' 

###################################################################### 

#  The  city  attractor  map  is  prepared  in  stages.  The  basic  approach  is 

#  as  follows: 

#  Fetch  the  small,  medium,  and  large  city  footprint  maps 

#  For  each  (small,  medium,  and  large) 

#  Create  road  travel  time  map  at  30  meter  resolution 

#  Create  overland  travel  time  map  starting  from  times  set  above 

#  Merge  the  two  maps 

#  Combine  the  small,  medium,  and  large  cumulative  maps 

###################################################################### 

#  Combine  the  travel-time  maps  to  get  to  nearest  small,  medium,  and  large 

#  city  into  one  map  using  an  inverse-distance  weighting  algorithm  using  the 

#  population  of  each  type  as  the  base  weight. 
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LEcities_att :  smallCityTime  mediumCityTime  largeCityTime 
@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' LEcities_att=eval (  \ 

lt=if (largeCityTime==0, 1 , largeCityTime)  /  6.0,  \ 
mt=if (mediumCityTime==0 , 1 , mediumCityTime)  /  6.0,  \ 
st=if ( smaliCityTime==0, 1 , smallCityTime)  /  6.0,  \ 
float  (${LARGEPOP) )/ (It  *  It)  +  \ 
float (${MEDIUMPOP}) / (mt  *  mt)  +  \ 
float (${SMALLPOP) )/ (st  *  st))' 

#  r.mapcalc  ' LEcities_att=round (LEcities_att) ' 

#  The  cumulative  cost  developed  at  30  meter  resolution  over  roads 

#  is  combined  with  the  grosser  90  meter  resolution  overland  times 
smallCityTime:  smallCityTime30  smallCityTime90 

@echo  '#########################' ;  echo  $@  ;  date 

g. region  res=30 

r.mapcalc  ' smallCityTime=\ 

if (isnull ( smallCityTime30 ) , smallCityTime90, smallCityTime30) ' 

#  The  cumulative  cost  developed  at  30  meter  resolution  over  roads 

#  is  combined  with  the  grosser  90  meter  resolution  overland  times 
mediumCityTime:  mediumCityTime30  mediumCityTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 

g. region  res=30 

r.mapcalc  ' mediumCityTime=\ 

if (isnull (mediumCityTime30) , mediumCityTime90, mediumCityTime30 ) ' 

#  The  cumulative  cost  developed  at  30  meter  resolution  over  roads 

#  is  combined  with  the  grosser  90  meter  resolution  overland  times 
largeCityTime:  largeCityTime30  largeCityTime90 

@echo  '  ######################### ' ;  echo  $@  ;  date 

g. region  res=30 

r.mapcalc  ' largeCityTime=\ 

if (isnull (largeCityTime30) , largeCityTime90, largeCityTime30) ' 

#  A  cumulative  cost  map  from  small  Cities  at  the  90  meter  resolution. 

#  The  results  averages  cell  values  across  90  meters  for  every  30  meter  cell 
smallCityTime90 :  smallCityTime30  travelTime90 

@echo  '  ######################### ' ;  echo  $@  ;  date 
r .mapcalc  ' tmp=if (isnull ( smallCityTime30) , 99999, small City Time 30 ) ' 
r.mapcalc  'tmp=min(\ 

tmp  [-1,-1],  trap  [-1,0],  tmp[-l,l],  \ 
tmp  [  0,-1],  tmp  ,  tmp [  0,1],  \ 
tmp  [  1,-1],  tmp  [  1,0],  tmp  [  1,1])  ' 

r.null  tmp  setnull=99999 
g. region  res=90 

r.cost  -v  -r  travelTime90  out=tmp2  start_rast=tmp  percent=100 

g. region  res=30 

r.mapcalc  ' smallCityTime90=  \ 

round ( (tmp2 [-1, -1 ]  +  tmp2[-l,0]  +  tmp2[-l,l]  +  \ 
tmp2  [  0,-1]  +  tmp2  +  tmp2 [  0,1]  +  \ 

tmp 2  [  1,-1]  +  tmp 2 [  1,0]  +  tmp2 [  1,1])  /  9)  ' 
g . remove  tmp, tmp2 

#  A  cumulative  cost  map  from  medium  Cities  at  the  90  meter  resolution. 

#  The  results  averages  cell  values  across  90  meters  for  every  30  meter  cell 
mediumCityTime90 :  mediumCityTime30  travelTime90 

@echo  ' ######################### ' ;  echo  $@  ;  date 

r .mapcalc  ' tmp=if (isnull (mediumCityTime30 ) ,  99  99  9, mediumCityTime30 )  ' 
r.mapcalc  'tmp=min(\ 

tmp [-1,-1],  tmp [-1,0],  tmp [-1,1] ,  \ 
tmp  [  0,-1],  tmp  ,  tmp [  0,1],  \ 
tmp  [  1,-1],  tmp  [  1,0],  tmp  [  1,1])  ' 

r.null  tmp  setnull=99999 
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g. region  res=90 

r.cost  -v  -r  travelTime90  out=tmp2  start_rast=tmp  percent=100 
g. region  res=30 

r.mapcalc  'mediumCityTime90=  \ 

round ( (tmp2 [-1, -1 ]  +  tmp2[-l,0]  +  tmp2[-l,l]  +  \ 
tmp2  [  0,-1]  +  tmp2  +  tmp2 [  0,1]  +  \ 

tmp2[  1,-1]  +  tmp2[  1,0]  +  tmp2 [  1,1])  /  9)  ' 
g . remove  tmp, tmp2 

#  A  cumulative  cost  map  from  large  Cities  at  the  90  meter  resolution. 

#  The  results  averages  cell  values  across  90  meters  for  every  30  meter  cell 
largeCityTime90 :  largeCityTime30  travelTime90 

@echo  1  ######################### 1 ;  echo  $@  ;  date 
r .mapcalc  ' tmp=if (isnull (largeCityTime30) , 99999, largeCityTime30 ) ' 
r.mapcalc  'tmp=min(\ 

tmp [-1,-1],  tmp [-1,0],  tmp [-1,1] ,  \ 
tmp  [  0,-1],  tmp  ,  tmp [  0,1],  \ 
tmp  [  1,-1],  tmp  [  1,0],  tmp  [  1,1])  ' 

r.null  tmp  setnull=99999 
g. region  res=90 

r.cost  -v  -r  travelTime90  out=tmp2  start_rast=tmp  percent=100 

g. region  res=30 

r.mapcalc  ' largeCityTime90=  \ 

round ( (tmp2 [-1, -1 ]  +  tmp2[-l,0]  +  tmp2[-l,l]  +  \ 
tmp2  [  0,-1]  +  tmp2  +  tmp2 [  0,1]  +  \ 

tmp 2  [  1,-1]  +  tmp 2 [  1,0]  +  tmp2 [  1,1])  /  9)  ' 
g . remove  tmp, tmp2 

#  A  cumulative  cost  map  at  the  30  meter  resolution  is  produced  out  to  the 

#  30-minute  travel  time  threshold.  With  a  sufficiently  fast  r.cost,  this 

#  step  could  suffice  for  the  entire  area. 
smallCityTime30 :  smallCity  travelTime30 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.cost  -v  travelTime30  out=smallCityTime30  start_rast=smallCity 
percent=100 

r.mapcalc  ' smallCityTime30=round ( smallCityTime30) ' 

#  A  cumulative  cost  map  at  the  30  meter  resolution  is  produced  out  to  the 

#  30-minute  travel  time  threshold.  With  a  sufficiently  fast  r.cost,  this 

#  step  could  suffice  for  the  entire  area. 
mediumCityTime30 :  mediumCity  travelTime30 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.cost  -v  travelTime30  out=mediumCityTime30  start_rast=mediumCity 
percent=100 

r.mapcalc  ' mediumCityTime30=round (mediumCityTime30) ' 

#  A  cumulative  cost  map  at  the  30  meter  resolution  is  produced  out  to  the 

#  30-minute  travel  time  threshold.  With  a  sufficiently  fast  r.cost,  this 

#  step  could  suffice  for  the  entire  area. 
largeCityTime30 :  largeCity  travelTime30 

@echo  1 ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.cost  -v  travelTime30  out=largeCityTime30  start_rast=largeCity 
percent=100 

r.mapcalc  ' largeCityTime30=round (largeCityTime30) ' 

#  Fetches  a  map  of  the  locations  of  the  small  cities 
smallCity : 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' smallCity=$ { SMALL_CITY } ' 
r.null  smallCity  setnull=0 
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#  Fetches  a  map  of  the  locations  of  the  medium  cities 
mediumCity : 

@echo  1  ######################### 1 ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' mediumCity=$ {MED_CITY } ' 
r.null  mediumCity  setnull=0 

#  Fetches  a  map  of  the  locations  of  the  large  cities 
largeCity : 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' largeCity=$ { LARGE_CITY } ' 
r.null  largeCity  setnull=0 

#  Fetches  a  map  of  no  growth  areas. 

#  Should  be  replaced  with  fetching  maps  of  various  no  growth  areas  that  can 
be 

#  edited  and  then  merged. 

LEno_growth : 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' LEno_growth=round ( $ { NO_GROWTH} ) ' 

#  Creates  a  slope  map  (in  degrees)  from  the  dem.  Assumes  the  dem  z  values 
are  meters 

LEslope:  dem 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r . slope . aspect  dem  slope=LEslope 
r.mapcalc  ' LEslope=round (LEslope) ' 

#  Fetches  the  digital  elevation  model 
dem: 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 
r.mapcalc  'dem=${DEM}' 

#  Create  a  travel  cost  surface  for  travel  to  ramps  (road  cat  6) 

LEramp_att:  travelTime30  travelTime90  road 

@echo  1 ######################### 1 ;  echo  $@  ;  date 

#  Start  the  travel  time  at  30  m  resolution 
g. region  res=30 

#  Look  for  ramp  (cat  6)  cells  that  border  limited  access  roads  (cat  1) 
r.mapcalc  tmp= ' if (road==6  &&  (\ 

road [-1,-1] ==1  M  road [-1,0] ==1  ||  road [-1, 1] ==1  ||  \ 

road [  0,-1] ==1  II  road [  0,1]==1  ||  \ 

road  [  1,-1]==}.  ||  road  [  1,0]  ==1  ||  road  [  1,  1]  ==1)  ,  1,  null  ( )  )  ' 

#  Run  r.cost  from  these  locations 

r.cost  -v  start_rast=tmp  input=travelTime30  output=tmp30  percent=100 

#  Run  travel  time  at  90  meter  resolution 
g. region  res=90 

r.cost  -v  -r  start_rast=tmp30  input=travelTime90  output=tmp90 
percent=100 

#  Combine 

g. region  res=30 

r .mapcalc  ' LEramp_att=round (if ( is null ( tmp30 ) , tmp90, tmp30 ) ) ' 

#  g. remove  tmp30,tmp90 

#  Create  a  travel  cost  surface  for  travel  to  state  roads  (road  cat  6) 
LEstaterd_att :  travelTime30  travelTime90  road 

@echo  '#########################';  echo  $@  ;  date 

#  Start  the  travel  time  at  30  m  resolution 
g. region  res=30 

#  Look  for  ramp  (cat  6)  cells  that  border  limited  access  roads  (cat  1) 
r.mapcalc  tmp= ' if (road==6 , 1 , null ( ) ) ' 
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#  Run  r.cost  from  these  locations 

r.cost  -v  start_rast=tmp  input=travelTime30  output=tmp30  percent=100 

#  Run  travel  time  at  90  meter  resolution 
g. region  res=90 

r.cost  -r  -v  start_rast=tmp30  input=travelTime90  output=tmp90 
percent=100 

#  Combine 

g. region  res=30 

r.mapcalc  ' LEstaterd_att=round ( if ( isnull ( tmp30 ) , tmp90 , tmp30 ) ) ' 
g. remove  tmp30,tmp90 

#  Create  a  travel  cost  surface  for  travel  to  intersections 
LEinter sect_att :  travelTime30  travelTime90  intersection 

@echo  '#########################' ;  echo  $@  ;  date 

#  Start  the  travel  time  at  30  m  resolution 
g. region  res=30 

#  Look  for  ramp  (cat  6)  cells  that  border  limited  access  roads  (cat  1) 
r .mapcalc  tmp= 'if (intersection, 1, nu 11 () ) ' 

#  Run  r.cost  from  these  locations 

r.cost  -v  start_rast=tmp  input=travelTime30  output=tmp30  percent=100 

#  Run  travel  time  at  90  meter  resolution 
g. region  res=90 

r.cost  -r  -v  start_rast=tmp30  input=travelTime90  output=tmp90 
percent=100 

#  Combine 

g. region  res=30 

r .mapcalc  ' LEintersect_att=round (if (isnull (tmp30) , tmp90, tmp30) ) ' 
g. remove  tmp30,tmp90 

#  Finds  state  (2)  and  county  (3)  road  intersections.  Create  a  map  of 

#  just  state  and  county  roads.  Expand  these  locations  by  one  cell  and 

#  then  thin  them  down  to  one  cell  (to  take  care  of  any  small  map  errors) . 

#  Then  find  cells  that  are  surrounded  by  more  than  two  cells  containing 

#  roads . 

intersection:  road 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' tmp=if ( ( road==2 )  | |  (road==3),  1,  0) ' 

r.thin  tmp  output=tmp 
r.null  tmp  null=0 

r.mapcalc  ' intersection=if (tmp==l , if (2<  \ 

( tmp [-1,-1]  +  tmp [-1,0]  t  tmp [-1,1]  +  \ 
tmp [0,-1]  t  tmp [  0,1]  t  \ 

tmp  [  1,-1]  t  tmp [  1,0]  t  tmp [  1,1]), 1,0))' 
g . remove  tmp 

#  Create  the  municipal  buffer  file  for  LEAM 
LEcity_buff:  municipalBoundary  LEboundary 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
r .mapcalc 

' LEcity_buf f=if (isnull (LEboundary) ,l,if(isnull (municipalBoundary) ,0,1)) ' 

#  Create  a  municipal  boundary  map  based  on  the  municipal  boundary  in  PERM 
municipalBoundary : 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r. buffer  in=$ {MUNICIPAL_BOUNDARY}  out=municipalBoundary  dist=1.5 
unit=miles 

#  Fetch  the  study  area  boundary  from  PERM 
LEboundary: 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r .mapcalc  ' LEboundary=if (isnull ( $ { STUDY_AREA} ) , null ( ) , 1 ) ' 
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#  Create  a  distance  surface  to  the  nearest  cell  containing  water. 
LEwater_att:  LElandcover  distCost 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r .mapcalc  tmp= ' if (LElandcover ==11 ,1,0) ' 

r.null  tmp  setnull=0 

r. buffer  tmp  output=LEwater_att 

distances=30,  60,  90,  12 0,150,  180,210,  240,  270,  300,  330,  360,  390,  42 0,450,  480,  51 0,5 
40, 570, 600 

r. mapcalc  LEwater_att= ' if ( isnull (LEwater^att) , 600 , (LEwater_att-l ) *30) ' 
g . remove  tmp 

#  Create  a  distance  surface  to  the  nearest  cell  containing  forest  - 
LEf orest_att :  LElandcover  distCost 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r. mapcalc  tmp= ' if (LElandcover==41  | |  LElandcover==42  | |  \ 

LElandcover==43  | |  LElandcover==91 , 1 , 0) ' 

r.null  tmp  setnull=0 

r. buffer  tmp  output=LEf orest_att 

distances=100,200,300, 400, 500, 600,700,800,900,100,1100,1200,1300,1400,1500,1 
600,1700,1800,1900,2000,2500,3000 

r .mapcalc  LEf ores t_att= 'if ( isnull (LEforest_att) ,3000, (LEf ores t_att- 
1)*100)' 

g . remove  tmp 

#  Computes  travel  time  to  the  nearest  road  -  assuming  an  off-road  travel 
speed 

#  of  .5  mph.  Uses  the  buffer  program  rather  than  the  cost  accumulation 
program . 

LEroad_att:  road 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 

r. mapcalc  tmp= ' if (road==4 , 1 , null ( ) ) ' 
r. buffer  tmp  output=LEroad_att 

distances=30, 60,90,120,150,180,210,240,270,300,330,360,390,420,450,480,510,5 
40, 570, 600, 630, 660, 690, 72 0,750, 780, 810, 840, 870, 900, 930, 960, 990 


#  Convert  distance  in  meters  to  time  in  minutes  -  assuming  5  mph 

#  60  min  mile  hr 

#  X  min  =  - *  -  * -  *  Y  meter 

#  hr  1609  meter  .5  mile 

# 

#  X  min  =  60  /  (1609  *  .5)  min/meter  *  Y  meter 

#  X  min  =  .0745  min/meter  *  Y  meter 

#  X  min  =  .0745  min/meter  *  Y  meter 

#  Max  is  1000  meter  =  75  min 

r .mapcalc  LEroad_att= 'if (isnull (LEroad_att ) ,  75 , round ( (LEroad_att- 


1)*30*.0745) ) ' 
g . remove  tmp 

#  Creates  a  travel  time  across  90  meter  cells  by  choosing  the  shortest  time 

#  from  among  the  9  associated  30  meter  cells.  It  preserves  travel  time 
along 

#  road  cat  1  (limited  access)  and  preserves  the  highway  buffer 
travelTime90 :  travelSpeed90 

@echo  ' ######################### ' ;  echo  $@  ;  date 
g. region  res=90 

#  90m/cell  *  lmile/1609meter  *  60min/hr  =  3 . 35538min . mile/cell . hr 

#  3 . 35538min . mile/cell . hr  /  x  miles/hr  =  Y  min/cell 

r .mapcalc  ' travelTime90=if (travelSpeed90>0, 3 . 35538/travelSpeed90, 
null () ) ' 

#  Computes  the  time  required  to  traverse  a  cell  in  the  N-S  or  E-W  axes 
travelTime30 :  travelSpeed30 
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@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

#  30m/cell  *  lmile/1609meter  *  60min/hr  =  1 . 1 18 4 6min . mile/cell . hr 

#  1 . 1 18 4 6min . mile/cell . hr  /  x  miles/hr  =  Y  min/cell 

r .mapcalc  ' t ravel Time 3 0=if ( travel Speed3 0>0 ,  1 . 11 84  6/ travel Speed30, 
null () ) ' 


#  Assigns  travel  speed  to  cells  based  on  road  and  landcover  types. 

#  Speeds  are  in  miles  per  hour  (mph) 

travelSpeed30 :  LEno_growth  LEboundary  hwyBuff  roadTravelSpeed30 
landTravelSpeed30 

@echo  '#########################';  echo  $@  ;  date 

g. region  res=30 

r. mapcalc  travelSpeed30= ' \ 

if ( isnull (LEboundary) ,  null(),  \ 

if (isnull (roadTravelSpeed30) , null  ( ) ,  \ 

if (roadTravelSpeed30 ! =0,  roadTravelSpeed30,  \ 

if (LEno_growth,  null(),  \ 

if (hwyBuff,  null(),  landTravelSpeed30) ) ) ) ) ' 


travelSpeed90 :  hwyBuff  roadTravelSpeed30  landTravelSpeed30 
travelableLandSpeed 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 
r .mapcalc 

tmp= 'if(isnull( landTravelSpeed30) , travelableLandSpeed, landTravelSpeed30 ) ' 
r. mapcalc  tmp='\ 
eval (  \ 

highRoad  =  max (  \ 

roadTravelSpeed30 [-1, -1 ] , \ 
roadTravelSpeed30 [-1, 0] , \ 
roadTravelSpeed30 [-1, 1] , \ 
roadTravelSpeed30 [  0,-1 ],\ 
roadTravelSpeed30  ,\ 
roadTravelSpeed30 [  0,1], \ 
roadTravelSpeed30 [  1,-1 ],\ 
roadTravelSpeed30 [  1,0], \ 
roadTravelSpeed30 [  1,1] ),\ 
highLand  =  max (  \ 
tmp [-1, -1] , \ 
tmp [-1,0] , \ 
tmp  [  — 1, 1]  ,  \ 
tmp [  0,-1 ],\ 
tmp  ,  \ 
tmp [  0,1] , \ 
tmp [  1,-1 ],\ 
tmp [  1 , 0 ] , \ 
tmp [  1,1]),  \ 
haveBuffer  =  if (  \ 

hwyBuff [-1,-1]==!  ||  hwyBuff [-1,0] ==1  I  I  hwyBuf f [ -1 , 1 ] ==1  ||  \ 
hwyBuff [  0,-l]==l  ||  hwyBuf f==l  I  I  hwyBuff[  0,1]==1  ||  \ 

hwyBuff [  1,-1]==1  ||  hwyBuff [  1,0]==1  ||  hwyBuf f [  1,1]==1,  1,  0),  \ 
if (highRoad  >  0,  highRoad,  \ 
if (haveBuffer  >  0,  null(),  \ 
if (highLand  >  0,  highLand,  null()))))  ' 


r. mapcalc  travelSpeed90= ' \ 
eval  ( \ 

haveBuffer  =  if  (  \ 

isnull (tmp [-1 , -1 ] )  ||  isnull ( tmp [ -1 , 0 ] )  ||  isnull (tmp [-1 , 1 ] )  ||  \ 

isnull (tmp [  0,-1])  ||  isnull (tmp)  ||  isnull (tmp [  0,1])  ||  \ 

isnull(tmp[  1,-1])  ||  isnull(tmp[  1,0])  ||  isnull(tmp[  1,1]),  1,  0),  \ 

if (tmp==70, 70, if (haveBuffer, null ( ) , tmp) ) ) ' 


roadTravelSpeed30 :  road 
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@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  roadTravelSpeed30= ' \ 
if (road==l, 70,  \ 
if ( road==2 , 50 ,  \ 
if (road==3, 40,  \ 
if (road==4, 30,  \ 
if ( road==5, 20 ,  \ 
if (road==6, 30,  \ 
if (road==7 ,10,0))))))) ' 

landTravelSpeed30 :  LElandcover 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  landTravelSpeed30= ' \ 
if (LElandcover==2 1 , 10 ,  \ 
if  (LElandcover==22, 10,  \ 
if (LElandcover==23, 10, null ())))'  \ 

travelableLandSpeed :  LElandcover 

r.mapcalc  travelableLandSpeed= '  \ 
if (LEno_growth ! =1  &&  \ 

LElandcover ! =11  &&  \ 

LElandcover ! =12  &&  \ 

LElandcover ! =91  &&  \ 

LElandcover ! =92  ,  1,  0)  ' 

#  Creates  a  cost  in  distance  map.  Every  30  meter  cell  is  30.  Used  by 

#  the  water  and  forest  attractor  map  processes 
distCost:  LEboundary 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  ' distCost=if ( isnull (LEboundary) ,  null(),  30)' 

#  Creates  a  map  of  1-cell  wide  buffers  along  limited  access  highways 

#  except  where  other  roads  intersect.  Problems: 

#  -  No  distinction  between  bridges  and  true  intersections 

#  -  Parallel  roads  in  the  next  cell  are  considered  access  points 
hwyBuff:  road 

@echo  '  ######################### ' ;  echo  $@  ;  date 
g. region  res=30 
r.mapcalc  hwyBuff='\ 

if (road  >  0,  0,  if ( road [-1 , -1 ] ==1  ||  road [-1, 0] ==1  ||  road [ -1 , 1 ] ==1  ||  \ 
road [  0,-1] ==1  ||  road [  0,1]==1  ||  \ 

road [  1,-1] ==1  ||  road [  1,0]==1  ||  road [  1,1] ==1,  1,  0))' 

#  Fetches  the  landcover  file  from  PERM  and  masks  it  with  the  study  area 
boundary 

LElandcover:  LEboundary 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 

r.mapcalc  LElandcover= ' if ( isnull (LEboundary) , null () ,  \ 
if (isnull ($ {LANDCOVER} ) , 0, $ {LANDCOVER} ) ) ' 

#  Fetches  the  road  file  from  PERM  and  masks  it  with  the  study  area  boundary 
road:  LEboundary 

@echo  '#########################';  echo  $@  ;  date 
g. region  res=30 
r .mapcalc 

road= 'if (isnull (LEboundary) , null ( ) , if ($ {ROADS }==0, null ( ) , $ { ROADS } ) ) ' 
r.null  road  null=0 

CITIES_ATTRACTOR_COM_IND : LEcities_att 

r.mapcalc  'CITIES  ATTRACTOR  COM  IND=graph (LEcities  att,  \ 
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0.0,  0,580,0,1191,0.0608,  2184,0.1514,  3656,0.2830,  5826,0.4010, 
9615,-6704,  15755,1.0,  68706, .9929,  118000,1.0)' 

CITIES_ATTRACTOR_OPENSPACE : LEcities_att 

r.mapcalc  ' CITIES_ATTRACTOR_OPENSPACE=graph (LEcities_att,  \ 

0,0.0,  2000,-0924,  4000,0.9322,  6000,1.0,  8000,0.8383,  10000,0.5656, 
20000,0.4888,  30000,0,40000,0.1448,  118000,0.0)' 

#CITIES_ATTRACTOR_RES :LEcities_att 

#  r.mapcalc  ' CITIES_ATTRACTOR_RES=graph (LEcities_att,  \ 

#  0,0.0,  580,0.0,  1191,0.0464,  2184,0.0880,  3656,0.0648,  5826,0.1637, 
9615,0.2693,  15755,1.0000,  68706,0.5567,  118000,0)' 

CITIES_ATTRACTOR_RES : LEcities_att 

r.mapcalc  ' CITIES_ATTRACTOR_RES=graph (LEcities_att,  \ 

0,0.0,  20000,  1.0000)  ' 

CT_RD_ATTRACTOR_RES : LEroad^att 

r.mapcalc  ' CT_RD_ATTRACTOR_RES=graph (LEroad_att,  \ 

0,1,200,0.4044,  400,0.2687,  600,0.2123,  800,0.3106,  1400,0.0, 
1600,0,1800,0,2000,0,2500,0,2677,0)' 

RAM  P_AT  T  RAC  T  0  R_RE  S : LEramp^att 

r.mapcalc  ' RAMP_ATTRACTOR_RES=graph (LEramp_att,  \ 

0,0.0339,  1,0.5525,  2,0.8295,  3,1,4,0.8033,  5,0.6142,  6,0.5015,  7,0.2427, 
8,0.3111,  9,0.2729,  10,0.2293,  11,0.2577,  12,0.1855,  13,0.1943,  14,0.1534, 
15,0.1310,  16,0.1261,  17,0.1455,  18,0.1411,  19,0.1641,  20,0.0872,  21,0.0361, 
22,  0.0467,  23,0.0,  695,  0)  ' 

RD_INTERSECTION_RES : LEinter sect^att 

r.mapcalc  ' RD_INTERSECTION_RES=graph (LEintersect_att,  \ 

0,1,1,0.4499,  2,0.3401,  3,0.2182,  4,0.1560,  5,0.1632,  6,0.1400,  7,0.0969, 
8,0.0797,  9,0.0,  679,0) ' 

T  RE  E_AT  T  RAC  T  0  R_RE  S : LEforest_att 

r.mapcalc  ' TREE_ATTRACTOR_RES=graph (LEf orest_att,  \ 

0,1,  200,0.1145,  400,0.0158,  600,0.0030,  800,0.0,  2769,0)' 

WATER_ATTRACTOR_RES : LEwater_att 

r.mapcalc  ' WATER__ATTRACTOR_RES=graph (LEwater_att,  \ 

0,0.8807,  200,1,400,0.9947,  600,0.6515,  800,0.7220,  1000,0.5147, 
1200,0.2106,  1400,0.0791,  1600,0.1166,  1800,0.0648,  2000,0.0251,  2500,0.0, 
8377,  0)  ' 


DEM_RES : LEslope 

r.mapcalc  ' DEM_RES=graph (LEslope,  \ 

0.00,  1.00,  5.00,  1.00,  10.0,  1.00,  15.0,  1.00,  20.0,  0.935,  25.0,  0.31, 
30.0,  0.195,  35.0,  0.16,  40.0,  0.105,  45.0,  0.03,  50.0,  0.00,  55.0,  0.00, 
60.0,  0.00,  65.0,  0.00,  70.0,  0.00,  75.0,  0.00,  80.0,  0.00,  85.0,  0.00, 

90.0,  0.00)  ' 

S  T_RD_AT  T  RACT  OR_RE  S : LEstaterd_att 

r.mapcalc  ' ST_RD_ATTRACTOR_RES=graph (LEstaterd_att,  \ 

0,1,875,0.7913,  1785,0.5048,  2744,0.4177,  3765,0.3408,  4871,0.3229, 
6099,0.5550,  7512,0.2844,  8000,0.0,  13060,0)' 

CT_RD_ATTRACT0RJ30M_IND : LEroad_att 

r.mapcalc  ' CT_RD_ATTRACTOR_COM_IND=graph (LEroad_att,  \ 

0,1,200,0.5967,  400,0.4657,  600,0.3524,  800,0.3173,  1000,0.2330, 
1200,0.0,  1400,0.0,  1600,0,1800,0,2000,0,2500,0,2677,0)' 

CT_RD_ATTRACTOR_OPENSPACE : LEroad^att 

r.mapcalc  ' CT_RD_ATTRACTOR_OPENSPACE=graph (LEroad_att,  \ 

0,1,200,0.4044,  400,0.2687,  600,0.2123,  800,0.3106,  1000,0.5916, 
1200,0.6238,  1400,0.1263,  1600,0,1800,0,2000,0,2500,0,2677,0)' 
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RAMP_ATTRACTOR_COM_IND : LEramp^att 

r.mapcalc  ' RAMP_ATTRACTOR_COM_IND=graph (LEramp_att,  \ 

0,0.9254,  1,1,2,0.7202,  3,0.4261,  4,0.3061,  5,0.1874,  6,0.1017,  7,0.0839, 
8,0.0714,  9,0.0666,  10,0.0781,  11,0.0627,  12,0.0277,  13,0.0200,  14,0.0207, 
15,0.0293,  16,0.0205,  17,0.0148,  18,0.0116,  19,0.0140,  20,0.0158,  21,0.0134, 
22,  0.0047,  23,0.00,  695,0)  ' 

RAMP_ATTRACTOR_OPENSPACE : LEramp_att 

r.mapcalc  ' RAMP_ATTRACTOR_OPENSPACE=graph (LEramp_att,  \ 

0.00,  0.33,  10.0,  0.49,  20.0,  0.45,  30.0,  0.43,  40.0,  0.4,  50.0,  0.39, 
60.0,  0.38,  70.0,  0.35,  80.0,  0.325,  90.0,  0.295,  100,0.14)' 

RD_INTERSECTION_COM_IND : LEinter sect_att 

r.mapcalc  ' RD_INTERSECTION_COM_IND=graph (LEintersect_att,  \ 

0,1,1,0.4236,  2,0.2289,  3,0.1423,  4,0.0838,  5,0.05377  6,0.0409,  7,0.0403, 
8,0.0259,  9,0.0,  679,0) ' 

RD_INTERSECTION_OPENSPACE : LEinter sect_att 

r.mapcalc  ' RD_INTERSECTI01SM3PENSPACE=graph (LEintersect_att,  \ 
0,1,1,0.7619,  2,0.5447,  3,0.3471,  4,0.2191,  5,0.1479,  6,0.0723,  7,0.0368, 
8,0.0274,  9,0.0281,  10,0.0354,  11,0.0233,  12,0.0177,  679,0)' 

TREE_ATTRACTOR_COM_IND : LEforest_att 

r.mapcalc  ' TREE_ATTRACTOR_COM_IND=graph (LEforest_att,  \ 

0,1,200,0.1357,  400,0.0058,  600,0.0011,  800,0.0,  2769,0)' 

TREE_ATTRACTOR_OPENSPACE : LEforest^att 

r.mapcalc  ' TREE_ATTRACTOR_OPENSPACE=graph (LEf orest_att,  \ 

0.00,  0.405,  10.0,  0.385,  20.0,  0.375,  30.0,  0.365,  40.0,  0.35,  50.0, 
0.335,  60.0,  0.305,  70.0,  0.285,  80.0,  0.27,  90.0,  0.255,  100,0.06)' 

WATER_ATTRACTOR_COM_IND : LEwater_att 

r.mapcalc  ' WATER_ATTRACTOR_COM_IND=graph (LEwater_att,  \ 

0,1,200,0.5305,  400,0.3258,  600,0.2480,  800,0.1734,  1000,0.0954, 
1200,0.0706,  1400,0.0328,  1600,0.0,  8377,0)' 

WATER_ATTRACTOR_OPENSPACE : LEwater_att 

r.mapcalc  ' WATER_ATTRACTOR_OPENSPACE=graph (LEwater_att,  \ 

0.00,  0.455,  10.0,  0.44,  20.0,  0.435,  30.0,  0.43,  40.0,  0.415,  50.0, 
0.405,  60.0,  0.37,  70.0,  0.345,  80.0,  0.33,  90.0,  0.315,  100,0.235)' 


DEM_COM_IND : LEslope 

r.mapcalc  ' DEM_COM_IND=graph (LEslope,  \ 

0.00,  1.00,  5.00,  1.00,  10.0,  1.00,  15.0,  1.00,  20.0,  0.935,  25.0,  0.31, 
30.0,  0.195,  35.0,  0.16,  40.0,  0.105,  45.0,  0.03,  50.0,  0.00,  55.0,  0.00, 
60.0,  0.00,  65.0,  0.00,  70.0,  0.00,  75.0,  0.00,  80.0,  0.00,  85.0,  0.00, 

90.0,  0.00)  ' 

DEM_OPENSPACE : LEslope 

r.mapcalc  ' DEM_OPENSPACE=graph (LEslope,  \ 

0.00,  1.00,  5.00,  1.00,  10.0,  1.00,  15.0,  1.00,  20.0,  0.935,  25.0,  0.31, 
30.0,  0.195,  35.0,  0.16,  40.0,  0.105,  45.0,  0.03,  50.0,  0.00,  55.0,  0.00, 
60.0,  0.00,  65.0,  0.00,  70.0,  0.00,  75.0,  0.00,  80.0,  0.00,  85.0,  0.00, 

90.0,  0.00)  ' 

ST_RD_ATTRACT0RJ30M_IND : LEstaterd_att 

r.mapcalc  ' ST_RD_ATTRACTOR_COM_IND=graph (LEstaterd_att,  \ 

0,1,875,0.2825,  1785,0.1750,  2744,0.1578,  3765,0.0,  13060,0)' 

ST_RD_ATTRACTOR_OPENSPACE : LEstaterd_att 

r.mapcalc  ' ST_RD_ATTRACTOR_OPENSPACE=graph (LEstaterd_att,  \ 
0,1,875,0.4310,  1785,0.4017,  2744,0.1961,  3765,0.1695,  4871,0.1052, 
6099,0.1548,  7512,0.1381,  9356,0.0100,  13060,0)' 
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#  Evolve  Landscape 
futureLanduse :  ATT RAC T 0 R_RE S 

r .mapcalc 

' futureLanduse=if ( (ATTRACTOR_RES*ATTRACTOR_RES*ATTRACTOR_RES*ATTRACTOR_RES) > 
rand (0.0, 1.0) , if ( rand ( 0 . 0 , 1 . 0 ) > . 9 , 2 1 ,  null ( ) ),null() )  ' 

#  r. mapcalc  futureLanduse=if ' (ATTRACTOR^RES  *  ATTRACTOR_RES  *  ATTRACTOR_RES 
/  4  >\ 

#  rand (0 . 0, 1 . 0) , 2 1, null () ) ' 

check : 

d.what. rast 

ATTRACTOR_RES , CT_RD_ATTRACTOR_RES , DEM_RES, RAMP_ATTRACTOR_RES ,  ST_RD_ATTRACTOR 
_RES, TREE_ATTRACTOR_RES , WATER_ATTRACTOR_RES , CITIES_ATTRACTOR_RES, LEno_growth 
, LElandcover, LEboundary 

altfutureLanduse :  ATTRACTOR_RES 

#  x=$ (shell  r. stats  -c  ATTRACTOR^RES  |  sed  -e  's/-[A  ]*/*/'  -e  's/$$/+\\/' 
-e  '  s/A\* . */0/ '  I  be  )  ;  \ 

#  echo  $$x  ;  \ 

r. mapcalc  futureLanduse=if ' (ATTRACTOR_RES  *  ATTRACTOR_RES  >\ 
rand ( 0 . 0, 1 . 0) , 2 1, null () ) ' 

#  Removes  all  index  files 
cleanlndex : 

g. remove  rast=$ { Indexf iles } 

#  Removes  all  the  temporary  files  used  in  the  creation  of  the  LEAM  inputs 
clean : 

g. remove  rast=$ { Other files } 

#  Removes  all  the  temporary  and  LEAM  input  files 
veryclean:  clean 

g. remove  rast=$ { LEf iles } 
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